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Abstract. We consider the spectrum of the Fibonacci Hamiltonian for small 
values of the coupling constant. It is known that this set is a Cantor set of 
zero Lebesgue measure. Here we study the limit, as the value of the coupling 
constant approaches zero, of its thickness and its HausdorfT dimension. We 
prove that the thickness tends to infinity and, consequently, the Hausdorff di- 
mension of the spectrum tends to one. We also show that at small coupling, 
all gaps allowed by the gap labeling theorem are open and the length of every 
gap tends to zero linearly. Moreover, for sufficiently small coupling, the sum 
of the spectrum with itself is an interval. This last result provides a rigorous 
explanation of a phenomenon for the Fibonacci square lattice discovered nu- 
merically by Even-Dar Mandel and Lifshitz. Finally, we provide explicit upper 
and lower bounds for the solutions to the difference equation and use them to 
study the spectral measures and the transport exponents. 
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1. Introduction 



1.1. Background and Motivation. It is always exciting to obtain a new con- 
nection between two different areas of mathematics. Here we prove several new 
results concerning the spectral properties of the discrete Schrodinger operator with 
Fibonacci potential, the so-called Fibonacci Hamiltonian, using methods from the 
modern theory of dynamical systems (uniformly hyperbolic and normally hyper- 
bolic dynamics). 

The Fibonacci Hamiltonian is a central model in the study of electronic prop- 
erties of one-dimensional quasicrystals. It is given by the following bounded self- 
adjoint operator in ^^(Z), 

[Hv,uj'fp]{n) — ^{n + 1) + -^^{n — 1) + Vx[i-a,i){nc( + u mod l)^{n), 

where F > 0, a = and cj e T = M/Z. 

This operator family has been studied in many papers since the early 1980's and 
numerous fundamental results are known. Let us recall some of them and refer the 
reader to the survey articles |D00[ ID07a| IS95j for additional information. 

The spectrum is easily seen to be independent of w and may therefore be denoted 
by Ey. That is, <j{Hv,ui) = for every w G T. Indeed, this follows quickly from 
the minimality of the irrational rotation by a and strong operator convergence. 
It was shown by Siito that Ey has zero Lebesgue measure for every V > 0; see 
jS89j . Moreover, it is compact (since it is the spectrum of a bounded operator) 
and perfect (because the irrational rotation by a is ergodic). Thus, Ey is a zero- 
measure Cantor set. This result was recently strengthened by Cantat |Can| who 
showed that the Hausdorff dimension of Ey lies strictly between zero and one. 

Naturally, one is interested in fractal properties of Ey, such as its dimension, 
thickness, and denseness. While such a study is well-motivated from a purely 
mathematical perspective, we want to point out that there is significant additional 
interest in these quantities. In particular, it has recently been realized that the 
fractal dimension of the spectrum is intimately related with the long-time asymp- 
totics of the solution to the associated time-dependent Schrodinger equation, that 
is, idtcj) = Hv,cj^; see |DEGTj . 

Fractal properties of Ey are by now well understood for large values of V. Work 
of Casdagli [Casj and Siito |S87j shows that for V > 16, Ey is a dynamically 
defined Cantor set. It follows from this result that the Hausdorff dimension and 
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the upper and lower box counting dimension of all coincide; let us denote 
this common value by dimSy. Using this result, Damanik, Embree, Gorodetski, 
and Tcheremchantsev have shown upper and lower bounds for the dimension; see 
[DEGTj . A particular consequence of these bounds is the identification of the 
asymptotic behavior of the dimension as V tends to infinity: 

lim dimSy • \ogV ^ log(l + V2). 

The paper |DEGTj also discusses some of the implications for the dynamics of the 
Schrodinger equation; let us mention |DT071 IDT08| for further recent advances in 
this direction for the strongly coupled Fibonacci Hamiltonian. 

By contrast, hardly anything about Sy (beyond it having Hausdorff dimension 
strictly between zero and one) is known for small values of V. The largeness of V 
enters the proofs of the existing results in critical ways. Consequently, these proofs 
indeed break down once the largeness assumption is dropped. The purpose of this 
paper is to fill out this gap by completely different methods. 

We would like to emphasize that quantitative properties of regular Cantor sets 
such as thickness and denseness are widely used in dynamical systems (see |N791 
INTO] ■ |PTj . |Mej ) and they have found an application in number theory (see |AsOO[ 
lAsOli IAs02| . |Cus| . |Ha| . [HI]'), but to the best of our knowledge, these kinds of 
techniques have never been used before in the context of mathematical physics. 

1.2. Statement of the Main Results. In this subsection we describe our results 
for small coupling V. Clearly, as V approaches zero, Hy^^^ approaches the free 
Schrodinger operator 

[i?oV'](n) = ^(^ + 1) + i^in - 1), 

which is a well-studied object whose spectral properties are completely understood. 
In particular, the spectrum of Hq is given by the interval [—2,2]. It is natural to 
ask which spectral features of Hv,uj approach those of Hq. It follows from Siito's 
1989 result |S89| that the Lebesgue measure of the spectrum does not extend con- 
tinuously to the case V = 0. Given this situation, one would at least hope that the 
dimension of the spectrum is continuous at ^ = 0. 

It was shown by us in |DG09a| (and independently by Cantat [Can]) that Sy is a 
dynamically defined Cantor set for V > sufficiently small (i.e., the small coupling 
counterpart to Casdagli's result at large coupling). A consequence of this is the 
equality of Hausdorff dimension and upper and lower box counting dimensions of 
Sy in this coupling constant regime. Our first result shows that the dimension of 
the spectrum indeed extends continuously to = 0. 

Theorem 1.1. We have 

lim dim Sy = 1. 

More precisely, there are constants Ci, C2 > such that 

1 - CiF < dim Sy < 1 - Ca^ 

for V > sufficiently small. 

We get Theorem |1.1| as a consequence of a connection between the Hausdorff 
dimension of a Cantor set and its denseness and thickness, along with estimates for 
the latter quantities. Since these notions and connections may be less familiar to 



4 



DAVID DAMANIK AND ANTON GORODETSKI 



at least a part of our intended audience, let us recall the definitions and some of 
the main results; an excellent general reference in this context is |PT) . 

Let C C M be a Cantor set and denote by / its convex hull. Any connected 
component of I\C is called a gap of C. A presentation of C is given by an ordering 
U = {Un}n>i of the gaps of C. If m G C is a boundary point of a gap U of C, we 
denote by K the connected component of /\(C/i U C/2 U . . . U C/„) (with n chosen so 
that Un — U) that contains u and write 

With this notation, the thickness t(C) and the denseness 0{C) of C are given 

by 

t(C) = supinf t(C,Z^,u), 61(C) = inf sup t(C,W, it), 

and they are related to the Hausdorff dimension of C by the following inequalities 
(cf. [PTl Section 4.2]), 

— - — — - < dimn C < -, — - — — - . 

log(2+^) - -log(2 + ^) 

Due to these inequalities, Theorem |1.1| is a consequence of the following result: 

Theorem 1.2. We have 

lim t(T,v) — 00. 
More precisely, there are constants C3 , C4 > such that 

for V > sufficiently small. 

Bovier and Ghcz described in their 1995 paper |BG| the then-state of the art 
concerning mathematically rigorous results for Schrodinger operators in i'^ (Z) with 
potentials generated by primitive substitutions. The Fibonacci Hamiltonian be- 
longs to this class; more precisely, it is in many ways the most important example 
within this class of models. One of the most spectacular discoveries is that, in this 
class of models, the spectrum jumps from being an interval for coupling = to 
being a zero-measure Cantor set for coupling V > 0. That is, as the potential is 
turned on, a dense set of gaps opens immediately (and the complement of these 
gaps has zero Lebesgue measure). It is natural to ask about the size of these gaps, 
which can in fact be parametrized by a canonical countable set of gap labels; see 
|BBG92| . These gap openings were studied in [B] for a Thue-Morse potential and in 
[BBG91| for period doubling potential. However, for the important Fibonacci case, 
the problem remained open. In fact, Bovier and Ghez write on p. 2321 of BG : It is 
a quite perplexing feature that even in the simplest case of all, the golden Fibonacci 
sequence, the opening of the gaps at small coupling is not knownf^ 

Our next result resolves this issue completely and shows that, in the Fibonacci 
case, all gaps open linearly: 

^There is a perturbative approach to this problem for a class of models that includes the 
Fibonacci Hamiltonian by Sire and Mosseri; see |SM89| and [OKl ISl89l ISM90I ISMSj for related 
work. While their work is non-rigorous, it gives quite convincing arguments in favor of linear 
gap opening; see especially |SM89I Section 5]. It would be interesting to make their approach 
mathematically rigorous. 
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Theorem 1.3. For V > sufficiently small, the boundary points of a gap in the 

spectrum depend smoothly on the coupling constant V. Moreover, given any 
one-parameter continuous family {Uv}v>o of gaps of Tiy , we have that 

V-s-O \V\ 

exists and belongs to (0,oo). 

Figure [T] shows a plot of the spectrum for small coupling: 




Figure 1. The set {{E, V) : E e^v, < V < 2}. 



The plot illustrates the results contained in Theorems |1.1 - 1.3 It also suggests 



that the limit limy_j.o depends on the chosen family of gaps. We have more to 
say about the value of the limit in Theorem [L6] below. It will turn out that its size 
is related to the label assigned to it by the gap labeling theorem. 
Our next result concerns the sum set 

Ey + Sy = {El + E2 '■ El , E2 Sy}. 

This set is equal to the spectrum of the so-called square Fibonacci Hamiltonian. 
Here, one considers the Schrodinger operator 

(2) 

[Hy ^p]{m, n) —^jj{m + 1, n) + ^(m — l.n) + if{rn, n + 1) + ■0("^, n — 

+ V {x[i-a,i){ma mod 1) + X[i_a,i)(nQ mod l))ip{m,n) 

in The theory of tensor products of Hilbert spaces and operators then 

implies that (7{Hy^) = Ey + Ey, see Section [ij This operator and its spectrum 
have been studied numerically and heuristically by Even-Dar Mandel and Lifshitz 
in a series of papers |EL06[ IEL07[ IEL08| (a similar model was studied by Sire in 
|Si89] ). Their study suggested that at small coupling, the spectrum of Sy + Sy is 
not a Cantor set; quite on the contrary, it has no gaps at all. 
Our next theorem confirms this observation: 

Theorem 1.4. For V > Q sufficiently small, we have that a{Wy') = Sy + Ey is 
an interval. 



Certainly, the same statement holds for the cubic Fibonacci Hamiltonian (i.e., the 
analogously defined Schrodinger operator in with spectrum Ey + Ey + Ey). 
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Notice that Theorem |1.4| is a consequence of Theorem |1.2| and the famous Gap 
Lemma, which was used by Newhouse to construct persistent tangencies and generic 
diffeomorphisms with an infinite number of attractors (the so-called "Newhouse 



phenomenon"), see Subsection 6.2 for details: 

Gap Lemma (Newhouse |N79[|N70j ). IfCi,C2 C are Cantor sets such that 

t(Ci)-t(C2) >1, 

then either one of these sets is contained entirely in a ga^ of the other set, or 

Ci n C2 7^ 0. 

Let us turn to the formulation of results involving the integrated density of 
states, which is a quantity of fundamental importance associated with an ergodic 
family of Schrodinger operators. We first recall the definition of the integrated 
density of states. Denote the restriction of Hv,uj to some finite interval A C Z with 
Dirichlet boundary conditions by Hy^. We denote by N{E, u, V, A) the number of 
eigenvalues of Hy^ that are less than or equal E. The integrated density of states 
is given by 

(1) N{E,V)^ lim -N{E,Lj,V,[l,n]). 

n— >-oo n 

We will comment on the existence of the limit and some of its basic properties 
in Section |4] One of the most important applications of the integrated density 
of states is the so-called gap labeling. That is, one can identify a canonical set 
of gap labels, that is only associated with the underlying dynamics (in this case, 
an irrational rotation of the circle or the shift-transformation on a substitution- 
generated subshift over two symbols), in such a way that the value of N{E, V) for 
E G M. \ Sy must belong to this canonical set. In the Fibonacci case, this set is 
well-known (see, e.g., |BBG921 Eq. (6.7)]) and the general gap labeling theorem 
specializes to the following statement: 

(2) {N{E, V): EeR\ Sy} C {{ma} : m e Z} U {1} 

for every V ^ 0. Here {ma} denotes the fractional part of ma, that is, {ma} — 
ma — \jna\ . Notice that the set of gap labels is indeed V^-independent and only 
depends on the value of a from the underlying circle rotation. Since a is irrational, 
the set of gap labels is dense. In general, a dense set of gap labels is indicative 
of a Cantor spectrum and hence a common (and attractive) stronger version of 
proving Cantor spectrum is to show that the operator "has all its gaps open." For 
example, the Ten Martini Problem for the almost Mathieu operator is to show 
Cantor spectrum, while the Dry Ten Martini Problem is to show that all labels 
correspond to gaps in the spectrum. The former problem has been completely 
solved [A J) , while the latter has not yet been completely settled. Indeed, it is in 
general a hard problem to show that all labels given by the gap labeling theorem 
correspond to gaps and there are only few results of this kind. Here wc show the 
stronger (or "dry") form of Cantor spectrum for the weakly coupled Fibonacci 
Hamiltonian and establish complete gap labeling: 



^For the purpose of this lemma, we also consider the two unbounded gaps in addition to the 
bounded gaps considered above, that is, the connected components of the complement of the 
convex hull of the Cantor set in question. 
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Theorem 1.5. There is Vo > such that for every V G (0, Vq], all gaps allowed by 
the gap labeling theorem are open. That is, 

(3) {NiE, V): E eR\ Sy} = {{ma} : m e Z} U {1}. 

Complete gap labeling for the strongly coupled Fibonacci Hamiltonian was shown 
by Raymond in }Raj . where he proves (|3| for V > 4. We conjecture that ^ holds 
for every V > 0. 

Let us return to the existence of the limit in Theorem |1.3[ As was pointed out 
there, the value of the limit will depend on the family of gaps chosen. Now that the 
gap labeling has been introduced, we can refine the statement. For to e Z \ {0}, 
denote by Um{V) the gap of Sy where the integrated density of states takes the 
value {ma}. 

Theorem 1.6. There is a finite constant C* such that for every m 'Z \ {0}, we 
have 

\U„,{V)\ Cm 

lim — -— — = r 

v^o \V\ \m\ 

for a suitable Cm G (0, C*). 

Our final set of results concerns the spectral measures and transport exponents 
associated with the operator family. We will give precise definitions and state- 
ments of our results in Section [5] and limit ourselves to a brief description here. 
The ultimate goal of any analysis of a given Schrodinger operator is always an un- 
derstanding of the associated unitary group, which then allows one to understand 
the dynamics of the associated time-dependent Schrodinger equation. The standard 
transport exponents capture the spreading of the quantum state in space. Most 
approaches to a study of these transport exponents proceed via (time-independent) 
spectral theory and link continuity properties of the spectral measure, associated 
to the initial state of the time evolution via the spectral theorem, to lower bounds 
for the transport exponents. In one space dimension, these continuity properties 
can in turn be investigated via an analysis of the solutions of the time-independent 
Schrodinger equation. Our goal is to carry this out for the weakly coupled Fi- 
bonacci Hamiltonian. Indeed, results of this kind are known, but the dependence 
of the quantities entering the estimates on the coupling constant had not been 
optimized. We revisit these approaches here and improve them to yield the best 
possible quantitative estimates at small coupling that can be obtained with current 
technology. Our results in Section [5] are likely not optimal and in particular do not 
approach the (known) zero-coupling values. We regard it as an interesting open 
problem to either prove or disprove that the dimension estimates for the spectral 
measures and the lower bounds for the transport exponents approach the values in 
the free (zero coupling) case as the coupling approaches zero. 

Some of the results of this paper were announced in |DG09b| . 

1.3. Overview of the Paper. Let us outline the remaining parts of this paper. 
In Section [2] we give some necessary background information and recall how the 
trace map arises in the context of the Fibonacci Hamiltonian and some of its basic 
properties. Moreover, since our aim is an understanding of weak coupling phe- 
nomena, we discuss the case of zero coupling. Section |3] is the heart of the paper. 
Here we regard the weak coupling scenario as a perturbation of zero coupling and 
study the dynamics of the trace map and consequences thereof for the structure 
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of the spectrum as a set. Section |4] considers the integrated density of states and 
proves complete gap labehng at weak couphng and our quantitative version of the 
hnear gap opening result. Spectral measures and transport exponents are studied 
by means of solution estimates in Section [5] Higher-dimensional models gener- 
ated by a product construction are discussed in Section [6j where we confirm some 
predictions of Even-Dar Mandel and Lifshitz. Since their model is based on the off- 
diagonal Fibonacci Hamiltonian and the main body of this paper (and most of the 
other mathematical works on the Fibonacci Hamiltonian) considers the diagonal 
Fibonacci Hamiltonian, we develop all the basic results for the off-diagonal model 
in Appendix |A] and explain there how our work indeed confirms the predictions for 
the original Even-Dar Mandel-Lifshitz product model. 

To assist the reader in locating the proofs of the theorems from the previous 



subsection, here is where they may be found: We prove Theorem 1.2 (which im- 



plies Theorem 1.1) in Subsection 3.4 Theore m |1.3| in Subsection 3.1 Theorem 1.4 
in Subsection |6.3[ Theorem 1.5 in Subsection |4.2[ and finally Theorem 1.6 in Sub- 
section 1431 



Acknowledgment. We are grateful to Mark Embree for generating the plots shown 
in Figures |4] and [5] 

2. Preliminaries 

2.1. Description of the Trace Map and Previous Results. The main tool 
that we are using here is the so-called trace map. It was originally introduced in 
[Kal iKKTj : further useful references include [BUjj |BR| iHMl iRo] . Let us quickly 
recall how it arises from the substitution invariance of the Fibonacci potential; see 
|S87) for detailed proofs of some of the statements below. 

The one step transfer matrices associated with the difference equation Hv,ujU — 
Eu are given by 

, , (E -Vx[i-a,i)iTna + uj laodl) -1 

Denote the Fibonacci numbers by {Ffc}, that is, i^o = -^i = 1 and Fk+i = Fk + Fk-i 
for fc > 1. Then, one can show that the matrices 

M^,{E)^(\ M,{E) 



and 

Mk{E) = TvfliFk,E) X • • • X Tv,o{l,E) for fc > 1 
obey the recursive relations 

Mk+i{E) ^ Mk-iiE)MkiE) 
for fc > 0. Passing to the variables 

Xk{E) = ^TrMfe(S), 

this in turn implies 

Xk+i{E) = 2xk{E)xk-i{E) - Xk-2{E). 
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These recursion relations exhibit a conserved quantity; namely, we have 

Xk+i{EY + Xk{Ef + Xk^i{Ef - 2xk+i{E)xk[E)xk-i{E) - 1 = — 
for every > 0. 

Given these observations, it is then convenient to introduce the trace map 
T : ^ M^ T{x, y, z) = {2xy - z, x, y). 
The following functioi:j^ 

G(x, y, z) — x^ + y"^ + z^ — 2xyz — 1 

is invariant under the action of T, and hence T preserves the family of cubic sur- 
face^ 

Sv = I (2;, y, z) & -.x^ +y'^ + z^ - 2xyz = 1 + ^ 

Plots of the surfaces 5*0.01 and 6*0.5 are given in Figures [2] and [3j respectively. 

It is of course natural to consider the restriction Ty of the trace map T to the 
invariant surface Sy- That is, Ty : Sy — >■ Sy, Ty — T\s^. Denote by fly the 
set of points in Sy whose full orbits under Ty are bounded. A priori the set of 
bounded orbits of Ty could be different from the non- wandering sel|^of Ty, but our 
construction of the Markov partition and our analysis of the behavior of Ty near 
singularities show that here these two sets do coincide. Notice that this is parallel 
to the construction of the symbolic coding in jCasj . 

Let us recall that an invariant closed set A of a diffeomorphism / : Af — Af 
is hyperbolic if there exists a splitting of the tangent space T^M = E^ © at 
every point a; e A such that this splitting is invariant under Df, the differential Df 
exponentially contracts vectors from the stable subspaces {E'^}, and the differential 
of the inverse, Df~^, exponentially contracts vectors from the unstable subspaces 
{i?"}. A hyperbolic set A of a diffeomorphism f : M ^ M is locally maximal if 
there exists a neighborhood U oi K such that 

A = fl /"({/). 

We want to recall the following central result. 

Theorem 2.1 ([Casj, |DG09a) . [Can]). For V ^ 0, the set ^y is a locally maximal 
hyperbolic set of Ty : Sy — >■ Sy. It is homeomorphic to a Cantor set. 

Denote by ty the line 

{[E-V E \ 

It is easy to check that iy C Sy. 

The second central result about the trace map we wish to recall was proven by 
Siito in [S87l. 



The function G{x, y, 2) is called the Fricke character, or sometimes the Fricke-Vogt invariant. 
^The surface So is called the Cayley cubic. 

point p £ M of a diffeomorphism / : A/ — > M is wandering if there exists a neighborhood 
0(p) C M such that f^{0) n O = for any k e Z\0. The non-wandering set of / is the set of 
points that are not wandering. 



10 



DAVID DAMANIK AND ANTON GORODETSKI 




Theorem 2.2 (Siito 1987). An energy E belongs to the spectrum of Hy^u if o-nd 
only if the positive semiorbit of the point { ^^^ , under iterates of the trace 

map T is bounded. 

In fact, as also shown by Siito in |S87| . the trace map can be used to generate 
canonical approximations of the spectrum, Ey. Namely, consider the following sets: 

4^>-{SeM: for (a;„,y„,z„) =r"(^,f,l), we have min{|a;„|, |y„|} < 1}. 
Then, we have ^ e(7^^^ ^ ^v, that is, 

n 4"^- 

riGZ+ 

Figure [i] shows a plot of the sets {{£,¥) : E € Il^yK < F < |} for values 
of n up to 20j^ These plots illustrate nicely both the linear gap opening and the 
fact that the size of a gap depends on its label; compare Theorems |1.3| and |1.6| 
To further document linear gap opening through numerics, Figure [5] zooms into a 
portion of Figure [4] near a point {E, 0) for an energy E where a gap opens; we have 
chosen E sa 0.7248. 

2.2. Properties of the Trace Map for V — 0. We will regard the case of small 
y as a small perturbation of the case ^ = 0. This subsection is devoted to the 
study of this "unperturbed case." 

Denote by § the part of the surface So inside the cube {|a:| < l,\y\ < l,|z| < 1}. 
The surface S is homeomorphic to S*^, invariant, smooth everywhere except at the 
four points Pi = (1, 1, 1), P2 = (-1, -1, 1), P3 - (1, -1, -1), and P4 - (-1, 1, -1), 
where S has conic singularities, and the trace map T restricted to S is a factor of 



This is also how the plot in Figure [ij was obtained. In fact, what is shown there is the set 

{iE,v) : E e T.^y°\ < y < 2}. 
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-T.5 -1 -:a a -25 1 1.5 2.E 



Figure 4. The sets {{E,V) : E e < F < f } through n = 20. 

the hyperbohc automorphism of T'^ = M'^/Z^ given by 

A{9,ip) = [9 + if, 9) (mod 1). 
The semiconjugacy is given by the map 

F : {9, (f) H> (cos 2n{0 + ip), cos 2tt9, cos 2tt(p). 

The map A is hyperbohc, and is given by the matrix ^ = ^| , which has 
eigenvalues 

fi= — - — and -/i = — - — . 

Let us denote by v", v" e the unstable and stable eigenvectors of A: 

Av" ^ /iv", A^^' = ||v"|| = ||v"|| = 1. 

Fix some small C > and define the stable (resp., unstable) cone fields on in 
the following way: 

(4) K; = {v e TpR^ : V = wV" + \v'\ > r>"|}, 

= {v G TpR^ : V = i;"v" + ^;''v^ \v'^\ > r>"|}. 
These cone fields are invariant: 
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Figure 5. The sets : E eJ:\^>} near (0.7248,0) for n = 15,20,25,30. 

Also, the iterates of the map A expand vectors from the unstable cones, and the 
iterates of the map expand vectors from the stable cones: 

Vveif" VneN |A"v| > — ^=^"|v|, 

y^eK' VneN |A-"v| > ^^=^"|v|. 
p yiTc^ 

The families of cones {K''} and {K^} invariant under A can be also considered on 

The differential of the semiconjugacy F sends these cone families to stable and 
unstable cone families on S\{Pi, P2, -P3, -P4}- Let us denote these images by {/C*} 
and {/C"}. 

Lemma 2.3 (Lemma 3.1 from |DG09a| ). The differential of the semiconjugacy DF 
induces a map of the unit bundle 0/ to the unit bundle of S\{Pi, P2, P3, P4}. 
The derivatives of the restrictions of this map to fibers are uniformly bounded. In 
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particular, the sizes of the cones in the families {/C*} and {/C"} are uniformly 
hounded away from zero. 

Finally, consider the Markov partition for the map : — > that is shown 
in Figure [6] (and which had already appeared in |Cas) : for more details on Markov 
partitions for two-dimensional hyperbolic maps see [PTi Appendix 2]). Its image 
under the map _F : — )■ § is a Markov partition for the pseudo-Anosov map 
T : § ^ §. 




Figure 6. The Markov partition for the map A. 



3. The Spectrum as a Set 

It is not hard to see that the line Iv is transversal to the stable manifolds of the 
hyperbolic set fiy for small values of V (see, e.g., |DG09al Lemma 5.5]). Therefore 
the intersection of Iv and W'{Vlv) (and, hence, Sy) is a dynamically defined 
Cantor set (see, for example, [T]). In this section we study the properties of this 
one-parameter family of Cantor sets. Namely, in Subsection |3.1| we prove that the 
size of a given gap in the Cantor set tends to zero linearly as the coupling constant 
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(the parameter) tends to zero. In Subsection 3.2 we use normally hyperbolic theory 
to introduce a normalizing coordinate system in a neighborhood of a singularity. 
Then in Subsection |3.3| the order of the gaps that is related to the dynamics is 
chosen. Roughly speaking, the longer it takes for a gap to leave the union of the 
elements of the Markov partition, the higher is the order of the gap. Next, in 
Subsection |3.5| this normalizing coordinate system is used to study the distortion 
properties of the transitions through a neighborhood of a singularity. Finally, in 
the last three subsections we bring all the pieces together and prove the distortion 



property that immediately implies Theorem 1.2 



3.1. Linear Gap Opening as the Potential is Turned On. Here we prove 
Theorem 11.31 

Consider the dynamics of T in a neighborhood of Pi — (1,1,1). Due to the 
symmetries of the trace map this will also provide information on the dynamics near 
the other singularities. Take tq > small and let OroiPi) be an ro-neighborhood 
of the point Pi — (1,1,1) in M'^. Let us consider the set Per2{T) of periodic points 
of T of period 2; compare Figures [7] and [8j 




I.CO 1Q5 




Figure 7. S'o.i 
and Per2{T) 
near (1, 1, 1). 



Figure 8. S'0.2 
and Per2{T) 
near (1, 1, 1). 



Lemma 3.1. We have 

Per2{T) = l^{x,y,z) : x G (-00, \) U {\,oo), y 

Proof. Direct calculation. 



2x- 1 



□ 



Notice that in a neighborhood Ui of Pi, the intersection / = Per2{T) n Ui is 
a smooth curve that is a normally hyperbolic with respect to T (see, e.g., Appen- 
dix 1 in |PT) for the formal definition of normal hyperbolicity) . Therefore, the 
local center-stable manifold Wi^^{I) and the local center- unstable manifold Wi^^{I) 
defined by 

Wt:c{I) = {P e Or,{Pi) : T"(p) e 0,0 (A) for aU n G N} , 
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WiZil) ^{pe Or„{Pi) ■■ T-'\V) e Or„{Pi) for all n e N} 
are smooth two-dimensional surfaces. Also, the local strong stable manifold 
Wi^^^{Pi) and the local strong unstable manifold W;™(Pi) of the fixed point Pi, 
defined by 

W^ZiPi) - {P e Wt:M) ■■ T^{P) ^ Pi as n ^ 00} , 
WZiPi) = {P e : ^ Pi as n ^ 00} , 

are smooth curves. 

The Markov partition for the pseudo-Anosov map T : S — S can be extended 
to a Markov partition for the map Ty : Sv Sv for small values of V. Namely, 
there are four singular points Pi = (1, 1, 1), P2 = (—1, —1, 1), P3 = (1, —1, —1), and 
P4 = (—1, 1, —1) of S. The point Pi is a fixed point of T, and the points P2, P3, P4 
form a periodic orbit of period 3. For small V, on the surface Sv near Pi there is a 
hyperbolic orbit of the map Ty — P|sv of period 2, and near the orbit {P2, P3, P4} 
there is a hyperbolic periodic orbit of period 6. Pieces of stable and unstable 
manifolds of these 8 periodic points form a Markov partition for Ty : Sy — > Sy- 
For V 0, the elements of this Markov partition are disjoint. Let us denote these 
six rectangles (the elements of the Markov partition) by Ry, Ry, . . . , Ry. Let us 
also denote Ry — U^^iRy. 

It is convenient now to consider Ty : Sy — ^ Sy since for Ty, each of the eight 
periodic points that were born from singularities becomes a fixed point. Due to the 
symmetries of the trace map, the dynamics of is the same in a neighborhood of 
each of the singularities Pi, P2, P3, and P4. Pick any bounded gap O in £y\W^{Qy). 
The boundary points of O belong to stable manifolds of two fixed points near one 
of the singularities Pi,P2,P3, or P4. Without loss of generality assume that those 
fixed points are in a neighborhood of Pi. 

The set of fixed points of in a neighborhood [/i of Pi is a smooth curve 
Fix(T^Or„(Pl)) = Per2{T) n Oro(-Pi); see Lemma O above. Each of the fixed 



points has one of the eigenvalues equal to 1, one greater than 1, and one smaller 
than 1 in absolute value. Therefore the curve Fix(r^, Oro(A)) is a normally hy- 
perbolic manifold, and its stable set M^^(Fix(r^, Oro(Pi))) is a smooth two di- 
mensional surface; see |HPS] . The strong stable manifolds form a C^-foliation of 
VF"(Fix(T6,0^„(Pi))); see [PSWl Theorem B]. 

If pv and (\y are two fixed points of Ty in Oro(Pi), then these points form 
the intersection of the curve Fix(r^, 0^0 (Pi)) with Sy and can be found from the 
system 

V ^ 2x-l 

Z = X 

+y'^ + z"^ - 2xyz = 1 + ^ . 
If we parameterize Fix(T^, O^o (-Pi)) as {x = t + l,z = t+ l,y = ^^}, we get 

(t + If + + (t + If - 2(< + l)2l±l = 1 + — , 

^ ' 2t + l^ ^ ^ ' 2t + l 4' 

or 

4<'* + + 91^+^ + 1 _ ^2 
{2t + lf ^ ^ T' 

Since for the function f(t) = , we have /(O) 1, /'(O) ^ 0, /"(O) > 

0, the distance between the points py and qy is of order \V\ for small values of V. 
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Lemma 3.2. Let W CM.^ be a smooth surface with a -foliation on it. Let ^ C W 
be a smooth curve transversal to the foliation. Fix a leaf L <Z W of the foliation, 
and denote P = L f] ^. Take a point Q ^ P, Q £ L, and a line £q C M.^ , Q G £o, 
tangent to W at Q, but not tangent to the leaf L. Suppose that a family of lines 
{^y}ve(o,Vb) given such that iy as V ^ 0, each line £v , V > 0, intersects 

W at two points py and qy , and py Q, qy ^ Q as V ^ 0. 

Denote by Lp^ and Lq^, the leaves of the foliation that contain py and qy , 
respectively. Denote py — Lp^, D ^ and qy = Lq^ H ^. Then there exists a finite 
non-zero limit 

AisX{py,qy) 

lim -r. — ^■ 

v-j-o dist(pv, qy) 

Proof. Since £o is not tangent to the leaf L, there exists a curve <ZW tangent to 
^0, transversal to the foliation, and such that Q E ^. 

Set p = Lp^, n = Lq^, n ^. Since the foliation is C^, there exists a finite 
non-zero limit 

(5) lim f^ig^'y; ^ 0. 

V^O dist(py,qy) 

Let us consider a plane 11 tangent to at Q, and let vr : 14/^ 11 be an orthogonal 
projection (well defined and smooth in a neighborhood of Q). It is clear that 

(6) lim d'^fiPv)M^v)) ^ , 

V^O dlSt(py,qy) 

Also, since £y Iq a,s V 0, we have 

(7) lim dist(^(py),^(gy)) ^ ^ 

v^o dist{py,qy) 

Finally, since tt sends the C^-foliation of to a C^-foliation on 11, the projection 
along this foliation from TT(£y) to £q is C^-close to isometry. In particular, 

(8) lim distKgv)^^(£v)) = 1. 
v^o dist(7r(py),7r(qy)) 



The statement of Lemma 3.2 follows now from □ 



Proof of Theorem \1.3\ Apply Lemma |3.2| to the surface W = 

W'{Fi^{T\OroiPl))). □ 

3.2. Choice of a Coordinate System in a Neighborhood of a Singular 

Point. Due to the smoothness of the invariant manifolds of the curve of periodic 
points of period two described in Section |3.1[ there exists a smooth change of 
coordinates $ : Oro(A) such that $(Pi) (0,0,0) and 

• <&(/) is a part of the line {x = 0, z — 0}; 

• ^i^focW) is a part of the plane {z = 0}; 

• ^{Wf£(I)) is a part of the plane {x = 0}; 

• ^{Wf£{Pi)) is a part of the line {y = 0, z = 0}; 

• HWociPi)) is a part of the line {x = 0,y = 0}. 
Denote / = $ o To 
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In this case, 

/A-i 0\ 

A = D/(0,0,0) =D($oro$-i)(0,0,0) = -10 , 

\ OA/ 

where A is the largest eigenvalue of the differential DT{Pi) : Tp^R^ — > TpiM^, 

DT{P,)^ri ], A=^^=M^ 
\0 1 / ^ 

Let us denote 6v = ^{Sv)- Then, away from (0, 0, 0), the family {6y} is a smooth 
family of surfaces, (3o is diffeomorphic to a cone, contains the lines {y = 0,z — Q\ 
and {x = Q,y — 0}, and at each non-zero point on these lines, it has a quadratic 
tangency with a horizontal or vertical plane. 

Due to the symmetries of the trace map, similar changes of coordinates exist in a 
neighborhood of each of the other singularities. Denote Or^ — Oro(A) U Or„{P2) U 

0.0(^^3) u 0,0(^4). 

Fix a small constant C > and introduce the following cone fields in M?: 

(9) Kl^iw^TpM?, v = v,j,+v, : |v,| > Cy^|v,j,|}, 

(10) K; = {vGTpM3, v-v,+v^, : |v,| > C-i|v,^|}, 

(11) = {v e TpR\ V = V, + vy, : |v,| > C^\\Yy,\}, 

(12) - {v e TpM^ V = V, +v,^ : |v,| > C^Vs^.l}- 

Lemma 3.3. There are ri € (0, tq) and mo G N such that the following holds. 

1. T'^°{io)riOri is a union of two connected curvesji and^2, and ^('ji),i = 1,2, 
is tangent to the cone field K"; 

2. F-i(r'"o(4 n S)) is tangent to the cone field i^" (defined by Q). 

Definition 3.4. We will call the rectangle (the element of the Markov partition) 
the opposite to singularities Pi and P3 , and we will call the rectangle the opposite 
to singularities P2 and P4 

Notice that T"^° {Iq S) consists of a curve that connects Pi with an unstable 
boundary of i?g, a finite number of curves that connect two unstable sides of an 
element of the Markov partition and such that F^^{£^i) is tangent to the cone field 
JsT", and a curve that connects some other singularity with an unstable boundary 
of its opposite rectangle. 

Now let us take r2 e (0,ri) so smah that [J^^_^T^{Or2) C O^- 

For small V, denote by Sy.o^j the bounded component of Sv\Or2- The family 
{Sy.Or }ye[o,Vo) '-'1 surfaces with boundary depends smoothly on the parameter and 
has uniformly bounded curvature. For small V, a projection iry '■ Sv,o,-2 ^ ^ is 
defined. The map ny is smooth, and if p € S, <? S §^,0^3 j ^^'^ ''^v{q) = then TpS 
and TqSv are close. Denote by ICy (resp., ICy) the image of the cone /C" (resp., 
/C") under the differential of tt^^. 

Denote Fy = tt^^ o F, Fy : P-^i^XOr^) Sy. Denote also *y = $ o 
Fy. Compactness and Mean Value Theorem type arguments imply the following 
statement. 
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Lemma 3.5. There are Vq > and C > such that the following holds. Suppose 
that a^b £ T'^\F~'^{0r2), ^ Taf^, vi, G Ti,T'^ . Then the following inequalities 
hold for all V e [0,Vb]- 

dist{Fv {a), Fv{b)) < Cdist{a,b), 

dist(a,6) < Cdist{Fv{a),Fv{b)), 

Z{DFv{va),DFv{vb)) < C{Z{va,Vb) + dist{a,b)), 

A^a,Vb) < C{Z{DFviva),DFv{vb))+dist{Fv{a),Fv{b))). 

Moreover, if^/y{a) and^>v{b) are defined, then 

dist(*v'(a),^'y(6)) < Cdist{a,b), 

dist(a,6) < C'dist(*y(a),*y(6)), 

Z{D-i'viva),D^v{vb)) < C{Z{va, Vb) + dist(a, bj), 

Z{va,Vb) < CiZ{D^v{va), D^vivb)) + dist{^via),^vib))). 

Finally, notice that if C and Vq are taken sufficiently small, then the cone fields 
on and K",K" respect each other in the following sense. Suppose that 
a e T^\F-^{Or.,) is such that 'J'y(a) is defined for V G [0,Vo], and Va e K^. 
Lemma |2 . 3| implies that 

On the other hand, if 6 G T'^XF^^^Or^) is such that 5'y(6) is defined and D^vM e 
then 

Vb G K^. 

3.3. Ordering of the Gaps. To estimate the thickness of a Cantor set from below 
(or the denseness from above), it is enough to consider one particular ordering of 
its gaps. Here we choose a convenient ordering of gaps in £v H ^{ftv) (which is 
affinc equivalent to Sy). 

The trace map Ty, V ^ Q, has two periodic points of period 2, denote them 
by Pi{V) and P[{V), and six periodic points of period 6, denote them by P2{V), 
P^iy)-, ^3(^). Pd.{y\ Pi(y), and Pi{V). In Section|3l]we showed that the distance 
between Pi{V) and Pl{V) is of order \V\. 

We can choose the notation (swapping the notation for Pi{V) and P[{V), and/or 
for P2{V) and f2(^) necessary) in such a way that the following lemma holds. 

Lemma 3.6. // V is small enough, the line £v contains points Bi{V) € n 
W'iPiiV)) and B2{V) e £v Ci H/^"(P2(1^)) such that every point of the line which 
is not between Bi{V) and B2{V) tends to infinity under iterates ofTy. 

Denote by ly the closed interval on ly between the points Bi{V) and B2{V). 
It is known that the set of points on ly with bounded positive semiorbits is a 
dynamically defined Cantor set; see |DG09a[ ICan| . We would like to estimate the 
thickness of this Cantor set. 

Lemma 3.7. There are mo G N, < Ci < C2 and Vq > such that for all 
V G [0, Vo] the following holds. 

1. T™°(ly) n On is « union of two connected curves 71 (V^) and 72(^), arid 
$(7i(V")), i = 1, 2, is tangent to the cone field K"; 
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2. F^7\^"o(lv)\OrJ IS tangent to the cone field K""; 

3. Ty°{\y) consists of a curve that connects Wi^^{Pi) with a stable boundary 
of Bp , a finite number of curves S,i{V), each of which connects two stable sides of 
an element of the Markov partition and is such that Fy (^i) is tangent to the cone 
field K"^ , a curve that connects Wi^^{Pi)^ i G {2,3,4}, with a stable boundary of 
an opposite rectangle to Pi, and some "gaps" between the curves mentioned above. 
The length of these "gaps" is between CiV and C2V for all small enough V. 



Proof. The statement holds for V — (see Lemma 3.3). For V positive but small 
enough properties 1. and 2. hold by continuity, and property 3. follows from the 
fact that the distances between finite pieces of strong stable manifolds of points py 
and qy is of order V, and these strong stable manifolds form the stable parts of 
boundary of the Markov partition for Ty. □ 



We will call the preimages (under T™") of the gaps defined in Lemma 3.7 the 
gaps of order 1. 

Definition 3.8. A smooth curve 7 C Sy is tangent to an unstable cone field if 
Fy^i'^\0r2) is tangent to , and $(70 0^) is tangent to K". 

Definition 3.9. A curve is of type one if it is tangent to the unstable cone field 
and connects opposite sides of unstable boundaries of some i?*. 

A curve is of type two if it is tangent to the unstable cone field and connects 
a point from Wf^^{Pi)^i € {1,2,3,4}, with an unstable boundary of an opposite 
element of the Markov partition. 



In this terminology. Lemma 3.7 claims that Ty°{lv) consists of two curves of 



type two, some curves of type one, and some gaps between them of size of order V. 

Lemma 3.10. An image of a curve of type one under is a union of a finite 
number of curves of type one, and of a finite number of gaps of length between CiV 
and C2V . An image of a curve of type two under is a union of a curve of type 
two, a finite number of curves of type one, and a finite number of gaps of length 
between CiV and C2V . 

Proof. The first part follows from the properties of the Markov partition and the 
fact that the distance between strong stable (strong unstable) manifolds that form 
the Markov partition is of order V, see Subsection |3.1[ An image of a curve tan- 
gent to an unstable cone field is a curve tangent to an unstable cone field. Also, 
T^i^fociPi)) ^ ^ioc(-f»)- Therefore the image of a curve of type two under is 
a curve which is close to a finite piece of a strong unstable manifold of Pi, so the 
second part follows. □ 

Suppose that the gaps of order k have already been defined. Consider the com- 
plement of all gaps of order not greater than fc on ly. It consists of a finite number 
of closed intervals. Let J be one of them. Consider the curve 
construction, it is either a curve of type one, or of type two. In either case, the 
image r^(T™"+^^''"^' ( J)) = T™"+^''(J) consists of some curves of type one or two, 
and some gaps of size ~ V. Let us say that the preimages of these gaps (under 
Ty'>+'^'') are gaps of order k + 1. It is clear that every gap in ly n Wifly) has 
some finite order. Therefore we have ordered all the gaps. 
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3.4. Distortion Property: Estimate of the Gap Sizes. Let us consider ly 
and some gap (Z ly of order n. A bridge that corresponds to this gap is a 
connected component of the complement of the union of all gaps of order < n next 
to the gap. There are two bridges that correspond to the chosen gap, take one 
of them, and denote it by 75. Denote also 7 = 70 U 75. Now let us consider 
Tg = Ty°+^^"-^\jG) and Tg = r"°+^'""^^(7s). By definition of the order n of 
the gap we know that 

C3V<\^< c,v 
|i b\ 

for some constants C3 and C4 independent of V. 

Proposition 3.11. There is a constant K > 1 independent of the choice of the 
gap and of V such that 

\1b\ \Tb\ |7s| 



Notice that Theorem |1 .2| immediately follows from Proposition |3. 11 



The rest of this section is devoted to the proof of Proposition 3.11[ which is 



completed in Subsection 3.7 



3.5. Dynamics Near Singularities. Here we prove several technical propositions 
on the properties of the trace map in the coordinate system constructed in Subsec- 
tion 3.2 The first two propositions are reformulations of |DG09al Proposition 1]. 
The first one claims that a certain unstable cone field is invariant. 

We will use the variables {x, y, z) for coordinates in R^. For a point p G M^, we 
will denote its coordinates by (xp, yp, Zp). 

Proposition 3.12. Given Ci > 0, C2 > 0, A > 1, there exists 6q = So{Ci,C2, X) 
such that for any d G (0, 60), the following holds. 
Let f :M.^ ^M.^ be a -diffeomorphism such that 

(i) ll/llc^ < Ci; 

(ii) The plane {z = 0} is invariant under iterates of f; 

(iii) \\Df{p) — A\\ < S for every p G'R^ , where 



'X 



-1 











1 














A 



is a constant matrix. 
Introduce the following cone field in M"^ ; 

(13) K; = {w^ T,M\ V = v,j, + V, : |v,| > C2^\\^^y\}. 

Then for any point p — (xp, yp, Zp), \zp\ < 1 we have 

Df{Kp)<ZK^^^y 



Notice that the choice of the cone field Kp here (in (13)) and below (in (18)) 
corresponds to the choice of the cone field in ([9]). 

The next proposition establishes expansion of vectors from the introduced un- 
stable cones under the differential of the map. 
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Proposition 3.13. Given Ci > 0, C2 > 0, A > 1, e G (Oj i); ?/ > there exists 
So = 5q{Ci,C2,\,£), No e n,N„ ^ No{Ci,C2,\,e,So), and C = C{r]) > such 
that for any S € (0, Sq), the following holds. 

Under the conditions of and with the notation from Proposition \ 3.12 suppose 
that for the points p — (xp, i/p, Zp) and q = (xq, z/q, Zq), the following holds: 

1. < Zp <l and Q < Zq < I; 

2. For some N > Nq both f^{p) and f^{q) have z-coordinates larger than 1, 
and both f^~^{p) and f^~^{q) have z-coordinates not greater than 1; 

3. There is a smooth curve 7 : [0, 1] — t- M.'^ such that 7(0) = p, 7(1) — q, and for 
each t e [0, 1] we have ^'{t) e ^^(t)' 

If N > Nq (i.e., if Zp is small enough), then 

(14) |i?/'^(v)|>A^(i-4-)|v| for any v G X^, 
and if Df^{v) = u = u^y + u^, then 

(15) \u,y\ < 2(5i/2|u,|. 
Moreover, «/ jv^j > ?7|vx2;|, then 

(16) |£'/''(v)| > CA5(i-*^Vl for each fc = 1, 2, . . . , iV. 

In particular, 

(17) length(/^(7)) > A^Ci-^e) iength(7). 

In order to establish the distortion property we need better control over the 
expansion rates. In Proposition 3.15 we improve the estimates given by (14 1 and 
(17 1. As a first step we show that, roughly speaking, if a point stays for iterates 
in a neighborhood where normalizing coordinates are defined then it must be A~^- 
close to the center-stable manifold of a curve of fixed points. 

Proposition 3.14. Given Ci > 0, C2 > 0, A > 1, there exist Sq — (5o(Ci, C2, A), 
No = iVo(Ci,C2,A,(5o) G and C** > C* > such that for any 6 G {0,Sq), the 
following holds. 

Let / : M"^ — )• M'^ be a -diffeomorphism such that 

(i) ll/llc^ < Cu 

(ii) The planes {z = 0} and {x — 0} are invariant under iterates of f ; 

(iii) Every point of the line {z — 0,x = 0} is a fixed point of f; 

(iv) At a point Q €z {z — 0, x = 0} we have 



DfiQ) 




(v) \\Df(p) — A\\ < S for every p e M.^ , where 

^A-i 0^ 
A = DfiQ) = ( 10 

A; 
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Introduce the following cone fields in ; 



(18) 


if" 

p 


= {v e Tpl 


i\ 


V 


= ^xy + : 


Wz\ 




Wxy\}, 


(19) 




= {v e Tpl 


^^ 


V 


= Vx+ ^yz ■ 


|v.| 


< O.OIA"^ 


Wyz\}, 


(20) 




= {v e Tpl 


I', 


V 


= 'Vx+ ^yz ■■ 


|v.| 


> C2J\Xp 


\Wyz\}, 


(21) 


J^cs 


= {v e Tpl 


l\ 


V 


= V2 + V^y : 


|v.| 


< O.OIA-^I 


|Vxj/|}- 



Suppose that for a finite orbit pa,pi,pz, ■ ■ ■ ^pn we have 



{Pi))x > 1, {Pl)x < 1, {Pn)z > 1, {PN-l)z < 1, 

and there are curves 70 and such that 70 connects po with W^^ (Q) and is tangent 
to both cone fields and K"^^, and 7jv connects pn with W'^'^{Q) and is tangent 
to both cone fields and K'^^ . 
Then 

C*A-^ < |(po).| < C"A-^, and 
C*A-^ < \{pnU < (7"A-^. 

Proof. Consider an orthogonal from pq to the plane {z = 0}, and denote its base by 
Po- There is a unique point Qq on the line {z = 0, a; = 0} such that p^ € VF**((5o)- 
Denote the line segment connecting po and p^ by ctq and set CTj = f^{(To), i = 
1,2,. ..,7V. 

Similarly, consider an orthogonal from pjv to the plane {a; = 0}, and denote its 
base by p*j^. There is a unique point Qjv on the line = 0, x = 0} such that 
S W^^{Qiq). Denote the line segment connecting p^^ and by pjv, and set 
ft = /-^+'(pjv), i = 0,l,2,...,7V-l. 

We have 

< |cro| < |(Ti| < . . . < |(Tjv-i| < |o-jv|, 1 < |o-jv| < A(l + S), 

< Ipivl < |pjv-i| < ... < < IpoI, 1 < IpoI < A(l + 5). 

Denote hk = 6ist{pk,Q)- Then we have 
(A-^ -min(5,Ci6fe))|<Tfe| < \(7k-i\ < {X'^ +mm{5,Cibk))Wk\, k = l,2,...,N, 
{X-^ -mm{5,Cibk))\pk\ < \pk+i\ < {X-^ + mm{6,Cibk))\pk\, fc = 0, 1, 2, . . . , iV - 1. 

Now we have 

bk = dist{pk,Q) 

< dist(Q, Oo) + \pk\+ dist(g, Qn) + IcTfel 

< |o-fe| + \pk\ + C3C2{\/\pM + \/W\), 

where C3 does not depend on iV. Indeed, the distance between and W**(Q) is 
bounded above by the length of the curve 70, and since 70 is tangent to the cone 
fields K" and K'^", its length is not greater than C2\/|o'o|. On the other hand, 
dist((5, Qo) is of the same order as that distance since the strong stable manifolds 
of fixed points form a C^-foliation of the plane {z = 0}. In the same way one gets 
an estimate dist{Q.Q?f) < C:iC2\/\pN\- 

Since we have the a priori estimates \ak\ < {I + 5){X — (5)"^+*^ and \pk\ < 
(1 + 6){X — S)~'^, we also have 

bk < Wkl + \Pk\ + 2^2(1 + S){X - 
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If fc < iV/2, then 

6fe < (A - Sy^il + (A - 5)-^+2fc + 2C(A - 5)-^/2+'=) 

< C"(A-(5)-^ 

If fc > N/2, then 

< (A - + (A - 5)^-2'^ + 2C(A - 5)^/2-'=) 
Therefore we have 

N 

kol < IctatI J](A-i+mm(<5,Ci6fc)) 

lN/2] N 

< A(i + 5)x-^ n (1 + - '5)"'') • n (1 + <^i^'(^ - '5)^/'"'') 

fc=l fc=[iV/2l + l 

< C**A^^. 

Also, 

N 

\<Jo\ > WN\l[{X-'-min{6,Cibk)) 
fe=i 

lN/2] N 

> x-^ n (1 ~ CiC'ix - sy'') ■ Yl {I CiC'ix - 5)^/2-'=) 

k=l fc=[Ar/2] + l 

> C*A-^. 

In the same way we get estimates for p^. □ 

Proposition 3.15. Given Ci > 0, C2 > 0, A > 1, there exist Sq = (5o(Ci, C2, A), 
No <eN,No = 7Vo(Ci, C2, A, So) e N, and (7 > such that for any S G (0, So), the 
following holds. 

Under the conditions of and with the notation from Proposition \ 3.14\ suppose 
that V e Tp„M3^ V e K^^^. Then \Df^^{^r)\ > CA^/2|v|. 



Proof. We will use the notation from Proposition 3.14 and its proof. 

Let us denote = Df''{v),k = 0,1,..., N, and Dk = \{vk)z\,dk = \{vk)xy\- 
Let us normalize v in such a way that do = 1. Since v € K^^ and |cro| > C*A~^, 
we have Do > C^X^^^"^, where C5 is independent of N. 

Denote 

/ TOl(p) ti{p)\ 

Df{p) = rn^ip) e{p) t^ip) . 

\sM S2{p) X{p)J 

We have 

mi{p) ti{p)\ /^t\ I'v{p)Vj; + mi{p)vy + ti{p)vz 
m2{p) e{p) t2{p) = ■m2{p)^x + e{p)^y + t2{p)^z 

Sl{p) S2{p) X{p) J \vj \Sl{p)v^ + S2{p)Vy + X{p)vz 

Since ||/||c2 < Ci, we also have \iy{p)\ < X^^ + Cidist{Q,p), \mi{p)\, |m2(p)|, 
1^2(^)1 < Cidist{Q,p), and |A(p)| > X — Cidist{Q,p). Furthermore, if p 
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belongs to the plane {z = 0}, then si{p) ~ S2ip) — 0. Therefore, for arbitrary p, 
we have |si(p)|, |s2(p)| ^ C'l^p- This implies that we have the following estimates: 



(22) 



dk+i < (1 + niin(^, Cibk))dk + min(^, C\bk)Dk 
Dk+i > {\ ~ mm{S,Cibk))Dk - min{S,Ci\ak\)dk 



Lemma 3.16. There exists k* such that > for all k < k* , and dk < Dk for 
all k > k* . 

Proof. Indeed, if Dk > dk, then 

dk+i < (1 + S)dk + SDk < (1 + 26)Dk 

and 

Dk+i > (A - 5)Dk - 5dk > (A - 25)Dk. 
Since A — 25 > 1 + 2(5, we have Dk+i > dk+i- □ 

We have the following preliminary estimates, li k < N/2, then bk < C"(A — (5)^*^; 
if fc > N/2, then bk < C"(A - (5)"^+^ Also, |crfe| < (A - 6)-'^+^ for each k = 
0, 1, . . . , A^. Notice that this implies that niLi(l + C'l^i) is bounded by a constant 
that is independent of N. And, finally, Dk > CsX-^^^iX - S)''; see |DG09a[ 
Lemma 6.1]. 

If Dk < dk (i.e., k <k*), then 

dk+i < (1 + Cibk)dk + CibkDk 
< (1 + 2Cibk)dk 

' k 

Y[{l + 2Cibk) do 



< 



<Ce, 

where Cg does not depend on k or N. Moreover, we have 
Dk+i > (A - Cibk)Dk - Ci\ak\dk 

> (A - Cibk)Dk - CiC'eiX - 5)-^+'= 

CiCe (A-<5)-^+'=' 



> (A - Cibk)Dk 1^1 - 

> (A - C,bk)Dk (^1 - 

> (A - Cibk)Dk (l - 



A - Cibk Dk 
C,Ce (A - 5)-^+^ 



X-Cibk C5A-^/2(A-<5)'= 
C5(A-C: 



> A'^+^Do 



> A'^+^Dn 



n(l-(C^iA-^)6,:) 

4=0 
fe 

n(l-(^^iA-^)6,:) 



.4=0 



1 - 



1 - 



CiC( 



C^iX-C: 



C^iX-C: 



> CjX'^+'Do, 

since for any C > and ^ G (0, 1), one has lim7v-).oo(l — C^^)^ = 1. 
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If dk < Dk (i.e., k > k*), then 

Dk+i > (A - Cibk)Dk - Ci\<Jk\dk 

> iX-Cibk-Ci\ak\)Dk 

> XDk{l - X-^Cibk - X-^Ci\<Jk\) 

k 

> X'^+'CjDo n (1 - A-i^ifefc - X-'Ci\<Jk\) 



> CsX'^+^Do, 

where Cg does not depend on N or k. 

Finally, \Df^^{v)\ > > CsX^ Do > CA^/2|v|. 



□ 



Below (in the proof of Proposition 3.18) we will also need an estimate on k* 
provided by Lemma 3.16 Namely, we claim that k* cannot be much larger than 
N/2. The formal statement is the following. 

Lemma 3.17. There is a constant Cg independent of N such that 

Proof. We know that Dk*-i < dk'-i < Cq. Therefore Cg > > C^X^' Dq > 

CyA'^* • C^X-^I^, so A'^-* < (CeCf iC5~i)A^/2. □ 

Now we are ready to formulate the statement that will be used to check the 
distortion property of the trace map. 

Proposition 3.18. Given Ci > 0, C2 > 0,6*3 > 0, A > 1, there exist 6q — 
(5o(C'i,C2,C3, A), iVo = A'^o(Ci, C2, C3, A, ^o) G I^; and C > such that for any 
S € {0,Sq) and any A > 0, the following holds. 

Under the conditions of and with the notation from Proposition 3.14 suppose 
that the curve 70 has a curvature bounded by C3. Suppose also that for the points 
p = {xp, yp, Zp) and q — {xq, yq, Zq), the following holds: 

1- p,q& 70; 

2. For some N > No both f^{p) and f^{q) have z-coordinates larger than 1, 
and both f^^^{p) and f^^^{q) have z-coordinates not greater than 1; 

3. dist(/^(p),/^(g)) = A. 

Denote pk = f'^ip), Qk = /"(?), fc = 0, . . . , TV. 

Let V e TpM.^ and w G TqM.^ be vectors tangent to 70. 

Denote v*. = Df''{v) and = Df''{w), k = 0,...,N. 
between and Wk . 

Then, 



Let ttfc be the angle 



(23) 



N N 

ak < CA, and } dist{pk,qk) < CA. 



fc=0 



fc=0 



Proof. First of all, notice that it is enough to prove Proposition 3.18 in the case 
when the points p and q are arbitrarily close to each other. Indeed, otherwise 
split the piece of the curve 70 between the points p and q into a large number of 
extremely small pieces. If for each of them the statement of Proposition |3.18| holds, 
then by subadditivity of the inequalities (231 it holds in general. 
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Denote by T the piece of the curve 70 between the points po and qo, and set 
Tfc = /'^(r), k = 0, 1, 2, . . . , iV. Denote /i^ = |rfc|. Due to the remark above we 
can assume that for any vector tangent to F, the value of k* is the same. From the 
proof of Proposition |3. 15] we see that for fc = 0, 1, . . . , fc*, we have /i^ < Cg/ioj a-nd 
A^fiN> > CA^/2Cg-Vfc-, so /ifc. < ic-'C6)^-^/^^iN < ChA-^/^A. 

On the other hand, if k > k* , then 

A « > > CyA^-'^Dfc- > CyA^-'^^Mfc, 

where we denote by the length of the projection of F^ to the z-axis (slightly 
abusing the notation). Therefore, < (2C.^^)A~^+'^A. 
It follows that we have 

N N 

fe=0 k=0 

k* N 

= '^tJ'k+ 

k=0 k=k* + l 

N 

<k*■Ce^lo+ J2 (2Cfi)A-^+'=A 

k=k*+l 

< k* ■C'^X-^^^CeA + Ci2A 

< CA. 

Notice that for any two vectors v,w G X™ 

/A-i 0\ 

Z(Av,ylv^f) < AZ(v,w), where A= 10, 

\ OA/ 

and if a linear operator B is ^-closc to A, then 

Z(Bv, Sw) < (A + 0(1 + e)Z(v, w) 

= (A + e(A + l + 0)^(v,w) 
< A(l + 2e)Z(v,w). 

Therefore we have 
and 

ak+i < X{l + 2Cibk)ak + Cuik, fc = 0,l,...,r. 

Since nfc=o(-'^ ~^ 2Ci6fc) < C13 for some C13 that is independent of k* and N, we 
have 

ttfc < (A'^ + A'^-i + . . . + A + 1) • (Ci3CiC6C3)mo < C14AV0, 
where C14 is also independent of k and N. In particular, 

afc. < C\iX'''pio < (Ci4C9)A^/2C-U-^/2^ < C15A 

and 

k' k* 

Y,ak<Y, ^^i4AVo < CieA^'Vo < C17A. 

fe=0 fc=0 
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Now denote 



K^- = {v e TpR^ : |v,| > 100A|v,^|}. 



p 

If V, w e ii'"'', then Z(Av, Aw) < A" V2z(v, w), and the same holds for any Hnear 
operator B which is i5-close to A. 

There exists m G N independent of N such that if for a vector v we have 
|v^| > \v^y\, then Df^{v) e 

Also 

afc.+„ = Z{Dnvk^),Dr{wk,)) < Ci5(A + 25)™A = dgA, 

and for k > k* + m, we have a^+i < X~^/'^ak + Cijj^k- 
Denote v = A-^/^ _ ;^-i/2_ 
If Ci/ifc < yoik, then 

Q!fe+i < A-^/^afe + z/afc = (A'^/^ ^ j.^^,^ ^ A-^/^afe. 
If Ci/ifc > z/afe, then 

afe+i < A-i/^afe + C/Xfc < (A-V^i/-^ + l)CiMfc. 
Since Ef=fe.+™Mfc < ELfe.+™(2Cf ^)A-^+'=A < C12A, this implies that 

^ afe < CisA 

k=k*+m 

and 

^Qjfe < ^ttfc + (afe.+i + . . . + afc->+m-i) + ^ a/e < CA, 

fe=0 /s=0 /c=fe*+m 

concluding the proof. □ 

3.6. Distortion Property: Preliminary Estimates. The main result of this 

subsection is the following statement: 

Lemma 3.19. There are constants i? > 0, Vb > 0, and k > such that for any 

V € (0, Vb) and N G N, the following holds. Suppose that 7 C {lv)\Ori is a 
connected curve of length not greater than k. Let the points p,q € ly be such that 
Ty{p) e 7 and Ty{q) G 7, and Vp and Vq be unit vectors tangent to 7 at points p 
and q. Then 

N 

iZ{DTi,{vp),DT^{vq))+dist{TUp),TUq))) < R- 

1=0 

Notice that i^^^(§\Or2) is a torus without small neighborhoods of the prcimages 
of the singularities, and we can define the following map 

(24) fy :F-1(§\0^J-^T2, fv=F-^o-KvoTvo-K-^oF = Fy'^oTvoFv. 
If V is small, Ty is C^-close to the linear automorphism A on its domain. 

Lemma 3.20. For Vq > small enough, there exists t £ (0, 1) such that for 

V e [0, Vb], p,q G T'^\F~^{Or2) "■''^d unit vectors Vp e K^, Vq e K^, we have 

Z{Dfv,pivp), Dfv,q{vq)) < tZ{vp, Vq) + 2||f y || ^2 dist(p, q). 
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Proof. If Vq is small, then Ty is C^-close to the linear automorphism A. In partic- 
ular, for any point p G T^\F^^{0r2) ^^"^ vectors Vi,V2 £ K^, 

Z{Dfv,pivi),Dfv,p{v2)) < tZ{vi,V2), 

where t E (0, 1) can be chosen uniformly for all V E [0, Vq] and p E T^\F^^(Or2)- 
Therefore we have 

Z{Dfv,pivp),Dfv,M) < Z{Dfv,p{vp),DfvM) + Z{DfvM,DfvM) 

< tZ{vp,Vg) + 2\\DfvM - DTvMW 

< tZ{vp,Vq) + 2||TV||c2dist(p, q), 

as claimed. □ 
Definition 3.21. For any points p,q and any vectors Vp,Vq define 

(25) 

Lemma 3.22. For p,q E T^\F^^{OrJ, p 7^ q, and vectors Vp E , Vq E K^, 
consider the function ^{p,q,Vp,Vg) defined by (25). Suppose that p and q belong to 
a curve that is tangent to the unstable cone field. Then 

d{Tv{p),fv{q), DTv^pivp), Dfv,q{vq)) < t^{p,q,Vp,Vq) + 2||TV||c2- 
In particular, if ^{p,q,Vp,Vq) > -^-j^l^, then 

^{Tv{p),Tv{q),Dfv,p{vp),Dfv,q{vq)) < ^^^S{p,q,Vp,Vq). 
Proof. We have 

ZiDfv^p{Vp),Dfv,q{Vq)) 



d (TV (p) , Tv (g) , DTv,p (wp) , DTv,q (wq) ) = 



< 



distiTv {p),Tv{q)) 
tZ{vp,Vq) + 2\\fv\\c2distip, q) 



dist(p, q) 
= td{p,q,Vp,Vq) + 2||TV||c2- 

If we also have ^{p, q, Vp, Vq) > ^^^1^, then 

td{P,q,Vp,Vq) + 2||fy||c2 < t^{p,q,Vp,Vq) + ^-^^ {p, q, Vp, Vq) 



l + t 

2 



d{p,q,Vp,Vq) 



Lemma |3.5| immediately implies the following statement. 



□ 



Lemma 3.23. Fix a small V >0. Suppose that a,b E T'^\F ^(Or2) ^''6 such that 
1'via) and *y(6) are defined, and Va E Taf^XF-^iOrJ, v^ E ThT^\F-'^{0 
Then, 

Ha, b, Va,Vb) < C2(;?(*v(a), ^v{b), D-^vM, D^v{vb)) + 1) 

and 

di^^via), ^vib),D'^v{va),D^v{vb)) < C^^ia, b, Va, + 1). 
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Since Ty is C^-close to the linear automorphism A on its domain for smah > 0, 
we can assume that for aU V E [0, Vo]: we have 1171/11,^2 < 10. 

Let C be the constant from Proposition 3.18 where C3 — ( 3^ + l)C^ was taken. 
Fix a small r > 0. Take n* e N such that 

(26) (^l±iy' C^(C + 1)<^, and CV"* > A^"*^^"^)- 

Now we are going to choose a neighborhood U of the set of singularities 
{Pi, P2, P3, P4} in such a way that if an orbit of a point leaves U, then it does 
not enter U for the next n* iterates. Also, we will choose a smaller neighborhood 
U* C U such that if a small curve is tangent to an unstable cone field (see Defini- 
tion 3.8) and intersects U*, then either it is entirely inside of U, or its iterates will 



continue to intersect U until they reach the opposite rectangle. Here is how we do 
that. 

For small rg G (0,r2), we denote 0°^ = 0^3, O^, = ^(0^3) n Or,, O^^ = 
T{01-^) n Or, for each i > 1, Oj.^ = T-\Ol.+^) n Or, for each i < 0, and 
U ~ U^_3^(0*^). We will take rg so small that the following property holds. 
If p e 5v is such that T^^P) ^ V but p ^ C/, then T"(p) ^ U for every n e N with 
n < 71* . 

For small e (0,r3), we denote = 0^4, = ^(O^J n O^,, Oj.^ = 
r(0;;i) n O^, for each i > 1, 0\^ ^ r-i(0;+^) n Or, for each i < 0, and 

00 

(27) t^* = U (0^4)- 

z— — 00 

We will take r4 so small that the following property holds. Suppose 7 is a curve 
on such that ^'1/(7) is defined, Fyij) C Or,, 7 is tangent to the unstable cone 
field, Fvij) n Or AC/ 7^ 0, and Fv{-f) n U* 7^ 0. Then there is fc G N such that 
Ty{Fv{j)) n C/ 7^ for aU natural n < k, and Ty(Fy(7)) intersects the opposite 
rectangle of the Markov partition. 

Lemma 3.24. There are Vq > and Ri > such that the following holds. Suppose 
that V is a non-zero vector tangent to the line ly at some point p E ly ■ Let iV G N 
be such that Tv{p) belongs to the bounded component of Sv\Or,- Then, 

N 

J2\\DTI.{v)\\<Ri\\DTP{v)\\. 

i=0 

Proof. If V is small enough, the vector DT'^°{v) is tangent to the unstable cone 
field. Let us split the orbit {r™o(p),T'"«+i(p), . . .,T^{p)} into several intervals 

{r"^'?), T™"+l(g), . . . , T'^^-\q)}, {T'^^ (q),..., T'^^-^q)}, {T^^ (q),..., T^(g)} 

in such a way that the following properties hold: 

(1) for each i — 1,2, . . . , s, the points T*'^^^{q) and T^'-{q) are outside of 0^2! 

(2) if {r'''(g),...,T'="+i-i(g)}nC/* 7^0, then {T*^- (g), . . . , r'='+i-i(g)} C U*; 

(3) for eac h i = 1, 2, . . . , s — 1, we have either fc^+i — ki > n* (where n* is chosen 
due to ([26|) or {T'''{q), . . . ,T''^+^~^q)} nU* 7^0. 

Such a splitting exists due to the choice oi U* C f/ in (27) above. 

Apply Proposition |3.18| t o those intervals in the splitting that are contained in 
U* . The choice of n* in (|26|) above guarantees for the intervals that do not intersect 
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U* uniform expansion of the vector. The first and the last interval may have length 
greater than n* , and then we have uniform expansion that "kills" the distortion 
added by the change of coordinates, or smaller than n*, but then they do not add 
more than a constant to the sum. As a result, the required sum is bounded above 
by a geometrical progression. □ 

Lemma |3.24| implies the following statement. 

Lemma 3.25. There are constants Ri > 0, Vq > 0, and ki > such that for any 
V G (0, Vb) and N £ N, the following holds. Suppose that 7 C {ly)\Ori 
connected curve of length not greater than Hi . Let the points p, g € ly be such that 
(p) e 7 and (q) G 7- Then, 

N 

^dist(Ti,(p),T^(g)) <i?i. 

1=0 

Finally, the choice of n* and Proposition |3.18| imply that the function 
d{Tl,{p)^^q),DTl,{vp),DT^{vq)) is uniformly bounded, and together with 
Lemma [3. 25| this proves Lemma [3.19| 



3.7. Proof of the Distortion Property. 



Proof of Proposition \3. 1 1\ Notice that we need to prove that 



log 



VirsllTGl 



is bounded by some constant independent of the choice of the gap and of V. There 
are points pa G 7g a-nd pb ^ 1b such that if vg is a unit vector tangent to the 
curve 7g at pQ, and vb is a unit vector tangent to the curve 73 at ps, then 



log 



irdlTBl 

\^b\\ig\ 



log 
log 

n+1 



\T^^\ig)\\ib\ 



nn+2 



Tr'ilB)\\lG\ 



\DT, 



■n+2 



ivG)\ \ 



\DT--^'{vb)\J 



^ (log|OTv|^r.(„^)(T^(PG))l -log|OTv^lz,TM.B)(^^(PB))l 

i=0 

(+1 

< J2 |l0g \DTv\DTUva)(TUPG))\ - log \DTv\DT^^,,){n{pB))\ 
i=0 

n+1 

< J2 \\DTv\DT^^,a){Tl,iPG))\ - \DTy\oT^^,^^inipB))\ 



We estimate each of the terms in this sum using 

Lemma 3.26. Suppose / : M" — > M" is a smooth map, a,b £ M", and Va G 

TqM",!;;, G TfaM" are unit vectors. Then, 

\\Df\M\ - \DfUb)\\ < \\f\\c^iZ{va,v,) + \a-b\). 
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Proof of Lemma \3.26[ 

\\Df\vAa)\ - \DfUb)\\ < \\DflM\ - \DfUm + \\DfUa)\ " \DfUb)\\ 
<\\Df{a)\\-\va~vi,\ + \\f\\c--\a-b\ 

< \\f\\c<K~Vb\ + \a-b\) 

< \\f\\c2{Z{Va,Vb) + \a-b\). 

□ 



Now, Proposition 3.11 follows from Lemma 3.19 □ 



4. The Integrated Density of States 

4.1. Definition and Basic Properties. Recall the definition of N{E,V) given 
in ^, 

N{E,V)= lim -N{E,uj,V,[l,n\). 

n— >-oo 71 

Proposition 4.1 (Hof, see jHol ) . For every {E,V) G M^, the limit in ([T]) exists 
for every a; G T and its value does not depend on it. 

The following proposition collects some well-known properties of the integrated 
density of states. 

Proposition 4.2. (a) The map 

RxR3 {E,V) N{E, V) G [0, 1] 

is continuous. 

(b) For every y G M, there is a Borel measure on called the density of states 
measure and denoted by dNy , such that 

N{E,V) = f X{-oo.E]dNv. 
Jr 

(c) The topological support of the measure dNy is equal to Sy . 

(d) The density of states measure is the uj-average of the spectral measure associated 
with Hv,uj and the vector Sq G ^^(Z). That is, for every y G M and every bounded 
measurable g :M. ^ M., 

(28) f gdNv^ f {So,giHv,Mdto. 

(e) We have 

E <-2 

N{E, 0) = <( i arccos (- f ) -2 < E < 2 

E>2. 

Proof, (a) This follows from (the proof of) Lemma 3.1 and Theorem 3.2 in |ASj . 

(b) For every F G M, the map 

M.3 E^ N{E,V) G [0,1] 

is continuous by (a) and non-decreasing by construction, and hence it is the distri- 
bution function of a Borel measure on M. 

(c) and (d) See |CFKSl Section 9.2]. 
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(e) This is folklore; see, for example, \LSl Theorem 1.1] and its discussion there for 
a simple derivation. □ 

4.2. Complete Gap Labeling. Here we prove the following result, which implies 
Theorem 1 1 . 5| since the transversality assumption holds for Vo > sufficiently small. 



Theorem 4.3. Suppose Vq > is such that for every V € (0, Vq] and every point 
in fly, its stable manifold intersects ly transversally. Then, for every V € (0, Vq], 
all gaps allowed by the gap labeling theorem are open. That is, 

{N{E, V): E eR\ Sy} = {{ka} : fc e Z} U {1}. 

Proof. Consider the preimages of the singularities of the trace map F^^{Pi),i = 
1,2,3,4, on the torus. They form a set of 4 periodic points of the hyperbolic 
automorphism ^ : — >■ T^, and the stable manifolds of those periodic points 
intersect the line {</) ~ 0} transversally at the points {ka ( mod 1 )}, {ka + ^ ( 
mod 1 )}, {fca + f ( mod 1 )}, and {ka + ^ + 1 ( mod 1 )}. 

The images of these points under the semiconjugacy F form the set of points on 
£o of the form (± cos(7rma), ± cos(7rmQ!), 1), and they correspond to the energies 
E G {± cos(7rmQ!), m G Z}. The integrated density of states for the free Laplacian, 
N{E, 0), takes the values {{ma} : to e Z} at these energies. 

After we increase the value of the coupling constant, each singularity splits into 
two periodic points, and each of the stable manifolds of the singularities splits into 
two strong stable manifolds of the periodic points. Every point between the stable 
manifolds has an unbounded positive semiorbit, and therefore the interval that 
those manifolds cut in the line £v corresponds to a gap in the spectrum. 

Due to the continuous dependence of N{E, V) on the coupling constant and the 
local constancy of N{-, V) in the complement of Sy, the integrated density of states 
takes the same value in the formed gap as at the energy that corresponds to the 
initial point of intersection of the stable manifold of singularity with Iq. □ 

4.3. More on the Asymptotic Gap Lengths. 



Proof of Theorem 1.6 Fix any to G Z\{0}. The integrated density of states of the 
free Laplacian takes the values {±ma} at the energies {±cos7rmck}. 

If TO = 2k, then these energies correspond to points with ^-coordinates ka{ 
mod 1) and ka + ^{ mod 1) on F~^{£q). 

If TO = 2fc+l, then these energies correspond to points with ^-coordinates ka+^{ 
mod 1) and ka + ^ + ^{ mod 1) on ^-^(4). 

Take one of these points, Qfc G F^^{£o)- Let P* be the singularity such that 
Qk e P-^iWiP*)). Denote T ^ ^-^(VFf "(P*)). Let M e N he the smallest 
number such that ^^*^(r) contains Qk- 

Then A'^'^ {F^^{£q)) intersects F at some point Z, and the distance from Z to 
the set of singularities is uniformly (in |to|) bounded from zero. 

Denote by P'{V) and P"{V) the periodic points on Sy near P* (of period 2 
or of period 6, depending on P*). Denote by F' = Fy\W{'{P'{V)) and F" = 

Fy\wriP"{v)). 



Then T^^{Fy\£v)) (recaU that Ty was defined by ^) intersects F' and F" 
transversally near the point Z and the curves F' and F" cut in Ty {Fy^{ly)) an 
interval ly whose length is of order V . In other words, limy_5.o exists and is 
uniformly bounded from zero and from above. 
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We also have 



where C{Qk) is bounded from zero and from above. 

Notice that \iA~'^'{T) intersects F-i(£o) at Qk, then ^-^■^(r) must be of length 
at least \k\a-'^. On the other hand, \A-'^{T)\ = a-^|r|. Therefore, a-^|r| > 
|A;|a~-^. Hence 

On the other hand, 

\Um{V)\ 

is bounded from above (since is bounded). Therefore we have 

^ - I'^o \v\ |f--(/,)| - '^^'^'^Wl' 

where C{Qk) is uniformly bounded from above. □ 

Remark 4.4. Notice that one can actually claim a bit more. Namely, for those 
gaps that appear away from the endpoints of the free spectrum, the corresponding 
constant Cm is uniformly bounded away from zero. This follows from the fact that 
away from the singularities, the differential of F has norm which is bounded away 
from zero. 

5. Spectral Measures and Transport Exponents 

5.1. Solution Estimates. In this subsection we study solutions to the difference 
equation and prove upper and lower bounds for them. Results of this kind are 

known (wc provide the references below) , but our purpose here is to obtain explicit 
quantitative estimates, as functions of the coupling constant, as they enter explicitly 
in the bounds on fractal dimensions of spectral measures and transport exponents 
and wc wish to prove the best bounds possible for the latter quantities to get an 
idea about their behavior at weak coupling. These applications will be discussed 
in the next two subsections. 

Denote the largest root of the polynomial — (2 + V)x — 1 by ay. For small 
F > 0, we have ay ~ + cV with a suitable constant c. Our goal is to prove 

the following pair of theorems (recall that (j{Hv^uj) — for every oj). 

Let M{n,m,co,E) be the standard transfer matrix, that is, the SL(2,M) matrix 
that maps (u(rn+l), u{m))^ to {u{n+l),u{n)y for every solution u of the difference 
equation Hv,ujU = Eu. 

Theorem 5.1. For every V > and every 

log[(5 + 2F)V2(3 + y)a^] 



(29) C > 



logv2±i 



2 

there is a constant C > such that for every u and every E e.T,v, we have 
(30) max \\M{n,m,iJ,E)\\ < CN^. 

0<\n-m\<N 
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Theorem 5.2. For every V > 0, 

log(l+(2+WF) , , log[(5 + 2^)1/2(3 + i/)a^] 



16. log ' log^ 



\ — ) 2 

there are constants Cg, C„ > such that for every u, L >\, and E S Sy, we have 
that every solution to the difference equation Hy^uU — Eu, which is normalized by 
|u(0)p + |u(l)p = 1, obeys the estimates 

(31) C.L-y' <\\u\\l<CuL^\ 

where the local P' norm |1 • ||l is defined by 

[i] 

ll"lli = EK")l' + (^- W)KW + i)P- 

n=l 

forL> 1. 



Theorem 5.1 is a quantitative version of |DL99a[ Theorem 3], which in turn was a 
generalization of [IT. Theorem 1]. Theorem 5.2 on the other hand, is a quantitative 
version of jPKLi Propositions 5.1 and 5.2]. The upper bound in Theorem 5.2 is a 



consequence of Theorem 5.1 while the lower bound in Theorem 5.2 will be extracted 
from the details of the proof of [DKL, Proposition 5.1]. 

Let us recall some notation from |DL99a| . For y > 0, G M, and a e {0, 1}, we 
write 

Tiy.E,a)^{^-^^ / 

For a word w = w\ . . . Wn G {0, 1}* , we then set 

M{V,E,w) ^T(y,E,w„) X ••• X T{V,E,wi). 

Next, define the words {s„}„>o inductively by 

So = 0, si = 1, Sn = s„_is„_2 for n>2. 

By the definition of these words, there is a unique u e {0, 1}^+ (the fixed point 
of the Fibonacci substitution) that has s„ as a prefix for every n > 1; namely 

u = 1011010110110 We denote the set of finite subwords of u by Wu, that 

is, Wu = {0, 1, 01, 10, 11, 010, Oil, 101, 110, . . .}. By uniform recurrence of u it suf- 
fices to consider the matrices M{V, E, w) for E £ Sy and w e yV„ when proving 
Theorems O and El 

With these quantities, we have 

M^ = MjV,E)^M{V,E,Sn) 

and consequently 

x„ = x^iV,E) = ^TTMiV,E,s„) 

for n > 1. 

Finding upper bounds for the norm of transfer matrices consists of three steps. 
The first step is to bound the special matrices M„, that is, M(F„,0,a; = 0,E). 
This is the objective of Lemma |5.3[ Then, in Lemma |5.5[ we use interpolation 
to bound the norm of the matrices M(n, 0,a; = 0,E). Finally, the case of general 
matrices M{n, m, uj, E) is treated using partition. This will then complete the proof 
of Theorem 15.11 
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We begin with the first step. Part (a) of the following lemma is due to Siito (see 
[S87[ Lemma 2]) and part (b) is a relative of [ITl Lemma 4.(ii)]: 

Lemma 5.3. (a) We have E £ Sy if and only if the sequence is hounded. 

Moreover, we have 

V 

\Xn\ < 1 + Y 

for every E £ and n > 1. 

(b) With some positive V -dependent constant C , we have 

\\M.n\\<Cal 

for every E G Sy and n>\. 

Proof. As pointed out above, part (a) is known. Let us prove part (b). By the 
Cayley-Hamilton Theorem, we have that — 2x„M„ + / = 0, that is. 

The recursion for the matrices Mn, Mn = M„_2M„_i, gives 

Putting these things together, we find 

M„ = M„_2M«-i = M„_2(22;„_i/ - A/-_\) = 2a;„_iAf„_2 - M'}^. 

Since we also have ||M„|j = HM^-*^!!, we obtain from this identity along with part (a) 
the estimate 

I1M„|| < (2 + y)||Af„_2|| + ||M„-3|| 

for every i? € Ey. 

Consider for comparison the recursion 

m„ = (2 + y)m„_2 + m„_3. 

It is clear that we have ||Af„|| < m„ if we consider the sequence {to„} generated 
by the recursion and the initial conditions mi — ||A/i||, TO2 — ||-^^2|1, = ll-^sl!- 
Any solution of the recursion is of the form 

m„ — Cix^i + C2X2 + C3X3 , 

where the xj are the roots of the characteristic polynomial x^ — {2 + V)x — 1. Thus, 
the definition of ay implies that the estimate in (b) holds. □ 

Remark. Since (^^^i)", we can infer that ||A'/„|| < CF,\+^ with (5 ^ as 

F — >■ 0. Since for V = 0, the norm of the transfer matrices grows linearly at the 
energies E = ±2, we cannot expect a better result in general. However, since the 
transfer matrices are bounded when V = and E g (—2, 2) it is reasonable to 
expect that there is a better bound for small values of V for most energies in the 
spectrum. 

Let us now turn to the second step, which is the interpolation of the estimates 
from the previous lemma to non-Fibonacci sites in the to = sequence. The follow- 
ing lemma is a relative of |IT( Lemma 5] . Even though the proof is closely related 
to that in [ITj . we give the details to clearly show where the improved estimate 
comes from. 
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Lemma 5.4. For n > 1 and k > 0, we may write 

with coefficients that also depend on n and E. Moreover, for every n > 1 and 
E € Sy, we have that 

4 

^IP^'^I < (5 + 2V^)(3 + y)L^-/2J. 

Proof. Consider the case fc = 0. Then, 

M„M„ = Ml = 2XnMn - / 
by Cayley-Hamihon and hence we can set 

p(l) _ p(2) _ p(3) _ r, p(4) _ _-, 

It follows that for i? e Sy, we have the estimate 

Consider now the case k = 1. Then, 

MnMn+l = Mn+2 = 2x„+iMn + Af„_i - 2a;„_i/ 
by |IT1 Lemma 6] and hence we can set 

PW=0, Pf =2x„+i, Pf) = l, Pf)=-2x„_,. 
For E E Sy, we therefore have 

4 

^\Pi'^\ <5 + 2V. 

Next, assume that the lemma holds for k G {!,...,£} and consider the case 
k = £ + 1. Then, using [IT! Lemma 6] one more time, we find 

= (p,(1\m„+,.i + p/!\m„+,_2 + P/!'iM„+,_3 + Pt\l) Mn+e 

— Pl-l^'^n+i+l + Pi^i{2Xn+l-2Mn+l — Mn+l-l) + 

P^^_{{2xn+e^iM„+e-i - I) + Pf\M^+, 

= Pi^iMn+l+l + P^j^iMn+l + P/^\M„+£_i + P^ ■*/, 



where we set 



p(i) _ p(i) 

p(2) _ 9„ p(2) , p(4) 

p(3) _ _p(2) , r,^ p(3) 
— (.— 1 ' l-i (^—1 

p(4) _ _p(3) 
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It follows that 

Ei^il\i<(3+^)Ei^Hi- 

By induction hypothesis, we obtain the desired estimate. □ 



Lemma 5.5. For every > and every ^ obeying (29), there is a constant C so 
that 

||M(n,0,w = 0,E)\\ < Cn'^ 

for every n > I, E Cz Sy. 



Before we give the proof of Lemma 5.5 let us recall some facts about the minimal 



representation of positive integers in terms of Fibonacci numbers; compare, for 
example, (HCBl iLel [Z]. Given n > 1, there is a unique representation 

K 

such that Uk G Z+ and Uk+i — nk > 2. This is the so-called Zeckendorf represen- 
tation of n (which, incidentally, was published by Lekkerkerker some twenty years 
before Zeckendorf). It is found by the greedy algorithm, that is, i^„^ is set to be the 
largest Fibonacci number less than or equal to n, n is replaced by n — , and the 
process is repeated until we have a zero remainder. The property Uk+i ~nk > 2 fol- 
lows from the way the algorithm works together with the recursion for the Fibonacci 
numbers. Uniqueness of this representation in turn follows from Uk+i — rifc > 2. 
By construction, we have 

(32) <n<F„^+i <2F„,,. 

We will also need a relation between n and K, that is, the number of terms in 
the Zeckendorf representation of n. The algorithm to compute K = K{n) is the 
following. Start with two words over the positive integers, 1 and 1. The next word 
is obtained by writing down the previous word, following by the word preceding 
the previous word, but this time with every symbol increased by one. Now iterate 
this procedure: 

1, 1, 12, 122, 12223, 12223233, 1222323323334,... 

Concatenate all the words to obtain a semi-infinite word over the positive integers: 

111212212223122232331222323323334. . . 

The n-th symbol in this semi-infinite sequence is K{n). Notice that the lengths of 
these words follow the Fibonacci sequence. Moreover, only every other word has 
an increase in the maximum value of its entries, relative to the previous one, and 
the increase is by one whenever it happens. This shows that 

(33) ^2(if-i) < " 

Proof of Lemm,a \5.5[ From the Zeckendorf representation oi n > 1, one finds that 

M {n, 0,uj = 0,E) = Mn^M„, ■ ■ ■ M„^ ; 

compare the beginning of the proof of |IT1 Theorem 1] . 
Let us prove by induction that with 
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we have for every fc > 1, 



(34) 



II Af„„M„, • • • M„J1 < C(5 + 2\^)'= V 



fe,-no+2(fe-l) 



{avbvY 



Consider first the case k = 1. From Lemmas |5.3| and |5.4[ we may infer tliat 



||M„„M„J| < |fW_„J . ||M„J| + |p(f_„J . ||M„,_i|| + |p(f_„J . ||M„,_2|| + \Pi'l 
< C(5 + 2F)(3 + F)L("i-"'')/2Ja^i 



ni-no I 



Thus, the estimate (34) holds for k = I 



Now assume that (34) holds for k G {!,...,£} and consider the case k = £ + 1. 



Using Lemma 5.4 again, we see that 



= II (M„„M„JM„, • • •M„,M„,^J| 
< |Pi'L„J-||M„,M„,...M„,M„,^J 



(3) 

ni- 

.(4) 



|M„,_2M„, •••M„,M„,^J| 



+ \Pn:-no\-\\Mn,---Mn,M^,^,\\ 



\£i -"1+2+2(^-1) 



ne+i 



< \ J2\Pnlno\j Ci5 + 2V)%-^+''iavby) 

< (5 + 2V)6"i^"«C(5 + 2VYby"'+^'^ {avbvT'+' 
= C{5 + 2t/)^+ife-"''+2((^+i)-i) {aybvr^' . 



We conclude that ^ holds for fc = £ + 1. 
Therefore, 



\\Min,0,uj = 0,E)\\ < C(5 + 2y)^6-""+'(^-^) (ay&y)' 



THE WEAKLY COUPLED FIBONACCI HAMILTONIAN 



39 



It follows that 



log||Af(n,0.w = 0.£;)|| 
limsup < Imisup 



log (C(5 + 2\^)^5-"°+'(^-'^ (avbvr" 



< limsup 



< limsup 



log n n^oo log n 

K log(5 + 2V) + 2K log bv + uk log (ay 6y ) 



logn 



, log(5 + 2V) + 2 ^ log b 



log n 



log(5 + 2^)1/^ + log by + log (av6v) 

log[(5 + 2V)i/^(3 + ^)ay] 
log ^ 



where we used ( 32 1 and ( 33 ) in the third step. Since these estimates are uniform 
in e Sy, the result follows. □ 

In the third step, we can now turn to the 



Proof of Theorem \5.1\ Suppose V > and C obeys ( 29 ) . We will prove the following 
estimate, 

(35) \\M{V,E,w)\\ < C\w\^ for every V > 0, E eT^v, and w G W„, 
from which (|30[) immediately follows. 



Let V,E,w as in (35) be given. As explained in jDL99al Proof of Lemma 5.2], 
we can write 

w = xyz, 

where y G {a, b}* is a word of length 2 and a;^ (the reversal of a;) and z are prefixes 
of u. Thus, 

(36) \\M{V,E,w)\\ < \\M{V.E,x)\\ ■ \\M{V,E,y)\\ ■ \\MiV, E, z)\\. 
It follows from |DL99a[ Lemma 5.1] that 

(37) \\MiV,E,x'')\\ = \\M{V,E,x)\\. 
Moreover, Lemma [53] yields 

(38) \\M{V,E,x^)\\<C\x\i 
and 

(39) \\M{V,E,z)\\ <C\z\'^. 



Since C > 1 and y has length 2, (|35[) follows from (36|-(39) 



□ 



follows from Theorem 5.1 The 



DKLj since it is not made explicit 



Proof of Theorem 5^ The upper bound in (31 
lower bound in (31 1 will be extracted here from 
there. 

Denote U{n) = {u{n + l),u{n))^ and consider the associated local ^^-norms 
||J7||i. Power-law bounds for \\U\\l correspond in a natural way to power-law 
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ri/2 



bounds for so let us discuss the former object. By considering squares ad- 

jacent to the origin and Gordon's two-block method, Damanik, Killip, and Lenz 
showed that 

see [DKLl Lemma 4.1]. We find 

lo^^^^^^ log (l + 

limmf ^ > 



logi^sn 16-log(^) 

uniformly in cj G T and E G Sv'- This, together with monotonicity, yields the 
asserted lower bound. □ 

5.2. Spectral Measures. Given > and ut E T, let us consider the operator 
Since this operator is self-adjoint, the spectral theorem associates with each 
-0 e £^{Z) a Borel measure that is supported by Sy and obeys 

{^,giHv,^)^) = f g{E)dfi.^{E) 



for every bounded measurable function g on Ey. It follows from Theorem 5.2 
that each of these spectral measures gives zero weight to sets of small Hasudorff 
dimension. More precisely, we have the following result. 

Corollary 1. For every > 0, every lj, and every 



a < 



log (l + (S+W) + 16 • i"^) + + 2^)'/'(3 + V)av] ' 

we have that Ht(,{S) = for every ip £ £^(Z) and every Borel set S with h"{S) = 0. 



Proof. This is an immediate consequence of Theorem 5.2 and [DKLj Theorem 1]. 

□ 

5.3. Transport Exponents. Absolute continuity of spectral measures with re- 
spect to Hausdorff measures is important because it implies lower bounds for 
transport exponents via the Guarneri-Combes-Last Theorem. Let us recall this 
connection briefly. Given Hy^uj and ijj as above, consider the Schrodinger equation 
idtipit) = HY^^'ip{t) with initial condition ?A(0) — Tp- Then, the unique solution of 
this initial-value problem is given by = e-^*^^-"V- To measure the spreading 
of a wavepacket, one considers time-averaged moments of the position operator. 
This is of course mainly of interest when "0(0) = "0 is well-localized and in fact, in 
one usually considers the canonical initial state -0 = <^o- 
One is interested in the growth of {{\X\^)){T) as T cx) for p > 0, where 

and the time average of a t-dependcnt function / is given by either 
(40) {f){T)^^£ f{t)dt 
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(41) ^^^^^^^fi ^^''^"^/W^*- 

We will indicate in the results below which time-average is involved. However, for 
compactly supported (or fast-decaying) ip^ all the results mentioned below hold for 
both types of time-average. 

To measure the power-law growth of {{\X\p)){T), define 

iog((|x|p)(r) 

/?,b)=hmsup ^^^^^ 

and 



T-i-oo P log T 

The Guarneri-Combes-Last Theorem (see |La[ Theorem 6.1]) states that if /j,^ 
is not supported by a Borel set S with h°'{S) = 0, then /3*(p) > a for every p > 0, 



where the time-average is given by (40). 

Corollary 2. For every V>0,ujET,Oy^t/j£ ^^(Z), and p > 0, we have 



log(l+ (2^) +16 -log (^) +16 log[(5 + 2^)1/2(3 + l/)ay] 



where the time-average is given by (401 



Proof. Immediate from the Guarneri-Combes-Last Theorem and Corollary [T] □ 

We also have the following estimate, which holds for a special initial state but 
which is better when p is large: 

Corollary 3. For every > 0, a; G T, and p > 0, we have 

log^ 31og[(5 + 2F)i/2(3 + y)a^] 



± ^ & 2 



log ^ + log[(5 + 2y)i/2(3 + V)av] p (log ^ + log[(5 + 2y)i/2(3 + V)av]) 



where the time-average is given by (41). 



Proof. This is a consequence of Theorem 5.1 and results of Damanik, Siito, and 



Tcheremchantsev (see |DST[ Theorem 1] and compare also |DT081 Theorem 1]). □ 

It can be shown that for compactly supported (or fast-decaying) ip^ the transport 
exponents I3^{p) are non-decreasing functions of p taking values in [0,1]. Thus, 
in this case it is natural to consider the limits of these quantities as p | and 
p \ oo and denote them by /3^(0) and /3^(oo), respectively. The bounds in the 
two previous corollaries imply estimates for /3^(0) and P^{oo) as well. However, at 
small coupling, a better estimate for /3^(cxd) is obtained via a different route: 

Corollary 4. We have 

lim fff (oo) = lim lim [if (p) = 1, 



uniformly in a; G T. Here, the time-average is given by (41) 
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Proof. Damanik, Embree, Gorodetski, and Tcheremchantsev showed for the Fi- 
bonacci Hamiltonian that /3^(oo) > dinig(Ey) for every > and w e T; see 
[DEGTl Theorem 3]. Since the right-hand side converges to one asV^O, the resuh 
follows. □ 

Remarks, (a) As pointed out above, while the papers |DEGT[ IDSTI IDT03| work 



defined using ( 40 1 



with the time-average (41), the results above carry over to transport exponents 



(b) We find that /3^(oo), considered as a function of V, extends continuously to 
zero. While we would expect this also for other initial states tp, it does not follow 
from Corollary [2] It would be interesting to extend this continuity result to all 
fast-decaying initial states or to exhibit one for which it fails. 

6. Consequences for Higher-Dimensional Models 

6.1. Lattice Schrodinger Operators with Separable Potentials. Let d > 1 

be an integer and assume that for 1 < j < d, we have bounded maps Vj : Z — )• M. 
Consider the associated Schrodinger operators on ^^(Z), 

(42) [Hji^Kn) = ij{n + 1) + ijj{n - 1) + V,-(ri)V'(n). 
Furthermore, we let : Z^ ^ M be given by 

(43) V{n) = Viin,) + --- + VdM, 

where we express an element n oil/ as n = (rii, . . . , n^) with rij € Z. 

Finally, we introduce the Schrodinger operator on Pil/-) with potential V , that 

is, 

d 

(44) [Hil}]{n) = VX" + ej) + -^{n - e^)) + F(n)V'(n). 

Here, Cj denotes the element n of 1/^ that has rij — 1 and nt — for k ^ j. 



Potentials of the form ( 43 1 and Schrodinger operators of the form ( 44 ) are called 
separable. Let us first state some known results for separable Schrodinger operators. 

Proposition 6.1. (a) The spectrum of H is given by 

a{H)^a{Hi) + --- + a{Ha). 

(b) Given ipi, . . . ,ipd G ^^(2); denote by Hj the spectral measure corresponding to 
Hj and -ij^j . Furthermore, denote by the spectral measure corresponding to H and 
the element ip of£'^(Z'^) given by ^[n) — tpi{ni) ■ ■ -ipdin-d)- Then, 

M = /^i * • ■ ■ * Md- 

Proof. Recall the definition and properties of tensor products of Hilbert spaces 
and operators on these spaces; see, for example, [RSi Sections II. 4 and VIII. 10]. 
It follows from |RS| Theorem 11.10] that there is a unique unitary map U from 
£2(Z)(g)- • •(g)£2(Z) {d factors) to e'^{Z'^) so that for € ^^(Z), the elementary tensor 
• -(^ipd is mapped to the element ip olP{'L'^) given by ^{n) = -01 ("-i) • • • '4'd{nd). 
With this unitary map U, we have 

d 

U*HU = ^ Id • • • «) Id (g) i7j (g) Id (g) • • • (g) Id, 



THE WEAKLY COUPLED FIBONACCI HAMILTONIAN 



43 



with Hj being the j-th factor. Given this representation, part (a) now foUows from 
[RS[ Theorem VIII. 33] (see also the example on |RS[ p. 302]). Part (b) follows from 
the proof of |RSl Theorem VIII.33]; see also [BS] and [S]. □ 

6.2. A Consequence of the Newhouse Gap Lemma. We have seen above that 
the spectrum of a product model is given by the sum of the individual spectra. If 
these individual spectra are Cantor sets and we want to show that their sum is not 
a Cantor set, the following consequence of the Gap Lemma is useful. 

Lemma 6.2. Suppose C,K cM.^ are Cantor sets with t{C) ■ t(K) > 1. Assume 
also that the size of the largest gap of C is not greater than the diameter of K , and 
the size of the largest gap of K is not greater than the diameter of C . Then the 
sum C + K is a closed interval. 

Proof. Denote min C = Ci , max C — C2, min K = ki, max K — Let us prove 
that 

C + K =[ci + fci,C2 + fca]- 

The inclusion "C" is obvious, so let us prove the inclusion "D." Take an arbitrary 
point a; € [ci+fci, C2+fc2]. Then, x G C+X if and only if G C+K-x = C-{x-K). 
Therefore, 

xeC + K ^ Cn{x- K) ^d). 

Since t(C) • t{x — K) = t{C) ■ t{K) > 1, the Gap Lemma implies that a priori 
there are only four possibilities: 

(1) the intervals [ci,C2] and [x — k2,x — fci] are disjoint; 

(2) the set C is contained in a finite gap of the set {x — K); 

(3) the set {x — K) is contained in a finite gap of the set C; 

(4) Cn{x-K)^%. 

But the case (1) contradicts the assumption x G [ci + fci, C2 + fc2], and the cases (2) 
and (3) are impossible due to our assumption on the sizes of gaps and diameters of 
C and K . Therefore, we must have C n (x — K) ^ and hence x G C + if . □ 

6.3. The Square Fibonacci Hamiltonian. Let us now discuss the (diagonal ver- 
sion of the) model studied by Even-Dar Mandel and Lifshitz, namely the operator 

[Hy il)\{m, n) —il>{m + 1, n) + ipijn — 1, n) + ■(/'(to, rt + 1) + %}j{m, 7i — 

+ V {x[i~a,i){ma mod 1) + X[i-a, !)("•" mod l))V'(m-, n) 



in £ (Z ). By Proposition 6.1 we have 



(45) a(ii|,')) -Ey + 



Proof of TheoremVLA By (45), it suffices to show that Ey + Ey is an interval for 



y > sufficiently small. By Theorem 
V G (0, Vb). For such F's, Lemma [6^ 



1.2 there is Vb > such that r(Ey) > 1 for 



applies with C ~ K = Yjv and yields 
Ey + Ey = [min Ey + min Ey , max Ey + max Ey] , 
as desired. □ 
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Appendix A. The Off-Diagonal Fibonacci Hamiltonian 

The purpose of this appendix is to give complete proofs of the basic spectral prop- 
erties of the off-diagonal Fibonacci Hamiltonian. This operator has been considered 
in many physics papers and is the basic building block for the higher-dimensional 
product models studied by Even-Dar Mandel and Lifshitz. The mathematics liter- 
ature on the Fibonacci model provides an exhaustive study of the diagonal model, 
and it was always understood that "analogous results hold for the off-diagonal 
model." For the reader's convenience, we make these analogous results explicit 
here and hence provide in this paper a complete treatment of the spectrum of the 
Even-Dar Mandel-Lifshitz product model at weak coupling. 

A.l. Model and Results. Let a, b be two positive real numbers and consider the 
Fibonacci substitution, 

S{a) = ah, S{b) = a. 
This substitution rule extends to finite and one-sided infinite words by concate- 
nation. For example, S{aba) = abaab. Since S{a) begins with a, one obtains a 
one-sided infinite sequence that is invariant under S by iterating the substitution 
rule on a and taking a limit. Indeed, we have 

(46) S''{a) = S''-^{S{a)) = S'^'^ab) = S''-\a)S''-\b) = S''-\a)S''-^{a). 

In particular, S''{a) starts with S''~'^{a) and hence there is a imique one-sided 
infinite sequence u, the so-called Fibonacci substitution sequence, that starts with 
(a) for every k. The hull Cla,b is then obtained by considering all two-sided 
infinite sequences that locally look like u, 

^a,b = {w e {a, b}^ : every subword of w is a subword of u}. 

It is known that, conversely, every subword of u is a subword of every co £ ila,b- In 
this sense, u and all elements of the hull lo look exactly the same locally. 

We wish to single out a special element of fla,b- Notice that ba occurs in u and 
that 5'^ (a) = aba begins with a and S^{b) = ab ends with b. Thus, iterating S"^ on 
b\a, where | denotes the eventual origin, we obtain as a limit a two-sided infinite 
sequence which belongs to Q^.b and coincides with u to the right of the origin. This 
element of fla,b will be denoted by cOs- 

Each u) e Slafi generates a Jacobi matrix acting in 

With respect to the standard orthonormal basis {(5„}„£z of ^^(Z), where (5„ is one 
at n and vanishes otherwise, this operator has the matrix 

/■•• ••• ••• \ 






UJ-1 







U)-l 





Wo 




















Wi 






and it is clearly self- adjoint. 
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This family of operators, {H^}^(zq^ ^, is called the off-diagonal Fibonacci model. 
Of course, the structure of the Fibonacci sequence disappears when a = b. In this 
case, the hull consists of a single element, the constant two-sided infinite sequence 
taking the value a = b, and the spectrum and the spectral measures of the associated 
operator are well understood. For this reason, we will below always assume that 
a ^ b. Nevertheless, the limiting case, where we fix a, say, and let b tend to a is of 
definite interest. 

Our first result concerns general properties of the spectrum of H^. For C M, 
we denote by dimz/S* the Hausdorff dimension of S and by dimes' the box counting 
dimension of S (which is then implicitly claimed to exist). 

Theorem A.l. Suppose a,b > and a ^ b. Then, there exists a compact set 
C M such that a{H^) = f, for every uj G i^a.b, 

(i) T,a,b has zero Lebesgue measure. 

(ii) The Hausdorff dimension dim/f T,a,b is an analytic function of a and b. 
(in) < dim/f T,a,b < 1- 

More can be said about the spectrum when a and b are close to each other: 

2 I I 2 

Theorem A. 2. There exists eo > such that if a,b > 0, a ^ b, and 2ab < 1 + 
[in other words, if a and b are close enough), then 

(iv) The spectrum b is a Cantor set that depends continuously on a and b in 
the Hausdorff metric. 

(v) For every small 5 > and every E € '^a.b, we have 

diuiH {{E -6,E + 6)n ^a,b) = dims {{E -6,E + 6)n S„,b) 

= dim^f T^a.b dims ^a,b- 

(vi) Denote a — dhnH Tja^b, then the Hausdorff a-measure of Ha^b is positive 
and finite. 

(vii) We have that 

Y^a.b + Sa^b = [minEa.b + mini;Q^b,maxEQ^b -|- max 



Given these results, and especially Theorem A. 2 (vii), we can confirm rigorously 
the observation made by Even-Dar Mandel and Lifshitz in |EL06[ IEL07] that the 
square Fibonacci Hamiltonian (based on the off-diagonal one-dimensional model) 
has no gaps in its spectrum for sufficiently small coupling. 

Next, we turn to the spectral type of H^. 

Theorem A. 3. Suppose a,b > and a ^ b. Then, for every lo G ^a.b, has 
purely singular continuous spectrum. 

Throughout the rest of this Appendix we will only consider a, 6 > with a ^ b. 
Theorems I A. I and | A. 2 1 are proved in Subsection |A.2| and Theorem |A.3| is proved in 



Subsection IA.31 

A. 2. The Trace Map and its Relation to the Spectrum. The spectral prop- 
erties of are closely related to the behavior of the solutions to the difference 
equation 

(47) cj„+iM„+i + a;„u„_i = Eu^. 
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Denote 



Then u solves (47) (for every n e Z) if and only if U solves 
(48) U^^TUn,E)Un-i, 
(for every n G Z), where 



1 /E -I 



LO,, \UJ„ 



Note that detr„(n, £■) = 1. Iterating (|48|), we find 

[/„ - Af^(n,S){/o, 

where 

M^{n,E)=T^{n,E) x ■■■T^{l,E). 

With the Fibonacci numbers {Ffc}, generated hy = Fi — 1, F^^i = F^ + 
for fc > 1, we define 



Xk=Xk{E) = ^TTM^^{Fk,E). 



For example, we have 



a 

p^ 1/^^ -iVl^^ -IV/"^ ^ (E^-Ea^-E9 -E^ + b^^ 

and hence 

E E^-a^-b^ E^-2Ea^-Eb^ 

"^ = 2a' = 2ab ' ^ 2a^b ' 

Lemma A. 4. We have 

(50) Xk+i = 2xkXk^i - Xk-2 
for k > 2. Moreover, the quantity 

(51) Ik = xl^-^ + xl + xl_^ - 2xk+iXkXk-i - 1 

is independent of both k and E and it is given by 

(a^+6^_ 
4a262 

Proof. Since uig restricted to {n > 1} coincides with u and the prefixes s^ of u of 
length Fk obey s^+i = s^Sk-i for fc > 2 by construction (cf. (46)), the recursion 
(50) follows as in the diagonal case; compare [DOOi lD07a[ IS87| . This recursion in 
turn implies readily that (51) is fc- independent. 

In particular, the x^'s are again generated by the trace map 

T{x,y,z) = {2xy- z,x,y) 

and the preserved quantity is again 

I{x, y, z) ^ x'^ + y"^ + — 2xyz — 1. 



THE WEAKLY COUPLED FIBONACCI HAMILTONIAN 



47 



The only difference between the diagonal and the off-diagonal model can be found 
in the initial conditions. How are Xi, Xq^ x_i obtained? Observe that the trace 
map is invertible and hence we can apply its inverse twice to the already defined 
quantity {x3,X2,xi). We have 

T^^{x,y,z) = {y,z,2yz - x) 

and hence, using ( [iQ] ), 

(a;i,a::o,a;_i) = T^'^{x3,X2,xi) 

., : ^ iEa" - E9 E^-a^-b'' E' 

T- 



2a?h ' 2ah ' 2a 

. E^-a^-b^ E_ ^{E^ - a2 _ ^2)^ ^3 _ 2Ea^ - Eb'' 



T- 



2ab '2a Aa?b 2a^b 

E^^a^-b^ E E' 



2ab ' 2a 2b 



E E E^ E^ 
,2- 



2a' 26' iab 2ab 
E E a^ + r 



,2a' 26' 2a6 
It follows that 

I{xk+i,Xk,Xk-l) = I{xi,Xo, X-i) 

4a2 462 ^ 4a262 8a262 
_ (a^+6^)^ 

4a262 ^ 

for every fc > 0. □ 

It is of crucial importance for the spectral analysis that, as in the diagonal case, 
the invariant is energy-independent and strictly positive when a ^ b\ 

Lemma A. 5. The spectrum of is independent of u and may be denoted by 
T.a,b- With 

CTfc = {SeM:|xfc|<l}, 

we have 

(52) Ea,b = Pi cTfc U CTfc+i. 

fe>i 

Moreover, for every E € T,a b o-nd k > 2, 

and for E ^ Sa,6; l^^fel diverges super- exponentially. 

Proof. It is well known that the hull fla,b together with the standard shift transfor- 
mation is minimal. In particular, every w € ila,b may be approximated pointwise by 
a sequence of shifts of any other w £ fia.fc- The associated operators then converge 
strongly and we get ^(i?^) C (j{Hqj). Reversing the roles of uj and w, the first claim 
follows. 
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So let Sa b denote the common spectrum of the operators H^^, u e ^a.b- We 
have < max{2a, 26}. Thus, Y^a.b Q [- max{2a, 2&}, max{2a, 26}] =: la^b- For 

E € la.b, we have that at least one of xi, xq belongs to [—1, 1]. This observation 



allows us to use the exact same arguments Siito used to prove (52) for the diagonal 



model in [S87] . 

The only point where care needs to be taken is the claim that ak is the spectrum 
of the periodic Jacobi matrix obtained by repeating the values ujg takes on {1 < n < 
Fk} periodically on the off-diagonals. This, however, follows from the general theory 
of periodic Jacobi matrices, which relies on the diagonalization of the monodromy 
matrix (which is M^,^ {Fk , E) in this case) in order to obtain Floquet solutions and 
in particular discriminate between those energies that permit exponentially growing 
solutions and those that do not. This distinction works just as well here, but one 
needs to use that the aj„'s that enter in the U^s are uniformly bounded away from 
zero and infinity. 

Thus, after paying attention to this fact, we may now proceed along the lines 
of Siito. Let us describe the main steps of the argument. Since at least one of 
xi, Xq belongs to [—1, 1], we have a result analogous to |S87| Lemma 2] with the 
same proof as given there. Namely, the sequence {xk}k>o is unbounded if and 
only if there exists k such that \xk\ > 1 and |a;fc+i| > 1. Moreover, we then have 
\xk+i\ > c^' for some c > 1 and all / > 0. This shows 

Cfc U (T/c+l = [J t^k+l- 
l>0 

Using now the fact that the F^ periodic Jacobi matrices with spectrum converge 
strongly to H^^^ , we obtain 



k>ll>0 k>l k>l 



since the spectra ak and cTk+i are closed sets. Thus, we have one inclusion in (52 1. 

Next, suppose E £ r\k>i '^k^'^k+i If fc > 1 is such that \xk\ > 1, then |xfe_i| < 1 
and jxfc+il < 1. Since we have 



4+1 + 4 + 4-1 - 2xk+iXkXk-i - 1 = ^"".Jg^"/ - 1, 



4a262 



this implies 



Xk — Xk+lXk-l ± ( 1 — xl_^_i — xl_i + xl_^_ixl_i + ^^2^2 ' ^ 



and hence 

\xk\ < \xk+iXk-i\ + ( (1 - 4+i)(l - 4-l) + 



which, using < 1 and |a;fe+i| < 1 again, implies the estimate (53) for E G 

nfc>i '^k U (Jk+i- We will show in the next subsection that the boundedness of the 
sequence {xk}k>o implies that, for arbitrary uj € i^a,bj no solution of the difference 
equation ( [47| is square-summable at +00. Consequently, such E's belong to Sa_b[^ 



'^This follows by a standard argument: If E ^ ^a,b: then (Hij — E) ^ exists and hence 
(H^ — E)~^So is an £^(Z) vector that solves (|47| away from the origin. Choosing its values for 
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This shows the other inclusion in ( 52 ) and hence estabhshes it. Moreover, it follows 



that (53 1 holds for every E £ f, 



Finally, from the representation (52) of T,a,b and our observation above about 
unbounded sequences {xk}k>0: ^6 find that \xk\ diverges super-exponentially for 
E ^ Sa.6- This concludes the proof of the lemma. □ 

Lemma A. 6. For every E eR, there is j{E) > such that 

hm -log||Af^(n,i;;)|| =7(^), 

uniformly in u) £ ^a.b- 

Proof. This follows directly from the uniform subadditive ergodic theorem; compare 
|DL99bllIfo] . □ 

Lemma A. 7. The set Za.b '■= {E € M : ^{E) = 0} has zero Lehesgue measure. 

Proof. This is one of the central results of Kotani theory; see [Ko| and also jD07bj . 
Note that these papers only discuss the diagonal model. Kotani theory for Jacobi 
matrices is discussed in Carmona-Lacroix |CL| and the result needed can be deduced 
from what is presented there. For a recent reference that states a result sufficient 
for our purpose explicitly, see Remling |Re| . □ 



Lemma A. 8. We have , 



Proof. The inclusion T,a.b 3 Za.b holds by general principles. For example, one can 
construct Weyl sequences by truncation when ^{E) = 0. The inclusion T,a.b ^ Za.b 
can be proved in two ways. Either one uses the boundedness of Xk for energies 
E G Tja.b to prove exphcit polynomial upper bounds for ||M^(n, £')|| (as in |IT] 
for Lo = ujs 01 in [DL99b| for general lo £ ^a.b), or one combines the proof of the 
absence of decaying solutions at +oo for E £ Tia^b given in the next subsection 
with Osceledec's Theorem, which states that ^{E) > would imply the existence 
of an exponentially decaying solution at +oo. Here we use one more time that C/„ 
is comparable in norm to (m„,u„_i)"^. □ 



Proof of Theorems A.l\ and \A.S\ The existence of the uniform spectrum T,a.b was 



shown in Lemma A. 6 and the fact that 'Sa,b bas zero Lebesgue measure follows 



from Lemmas lA. 71 and lA. 81 The set of bounded orbits of the restriction of the trace 
map T : K3 £3 ^.j^g invariant surface I{x,y,z) = C = ^^^^^ - 1,C > 0, 
is hyperbolic; see |Can| (and also ^DG09a] for C sufficiently small and [CasJ for C 



sufficiently large). Due to Lemma A. 5 the points of the spectrum correspond to 
the points of the intersection of the line of the initial conditions 

E_ E_ a^^h^ \ 
2a 26' 2ab J ' ^ 

with the stable manifolds of the hyperbolic set of bounded orbits. Properties (ii) and 
(iii) can be proved in exactly the same way as Theorem 6.5 in [Can) . The line £a.b 
intersects the stable lamination of the hyperbolic set transversally for sufficiently 
small C > 0, as can be shown in the same way as for the diagonal Fibonacci 
Hamiltonian with a small coupling constant; see |DG09a| . Therefore the spectrum 
Sa^fc for close enough a and 6 is a dynamically defined Cantor set, and the properties 



n > 1, say, and then using l |47| to extend it to all of Z, we obtain a solution that is square- 
summable at +oo. 
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(iv)-(vi) follow; see |DEGT1 IDG09a| iMal iMMl E TT\ and references therein. The 
statement (vii) follows as in the diagonal case since the thickness of b tends to 

2 I .2 

infinity as approaches 1. □ 

Notice that a proof of the transversality of the line £a,b to the stable lamination of 
the hyperbolic set of bounded orbits for arbitrary a ^ b would imply the properties 
(iv)-(vi) for these values of a and b. 

A. 3. Singular Continuous Spectrum. In this subsection we prove Theorem |A.3| 

Given the results from the previous subsection, we can follow the proofs from the 
diagonal case quite closely. 

Proof of Theorem \A.!^ Since the absence of absolutely continuous spectrum follows 
from zero measure spectrum, we only need to show the absence of point spectrum. 
It was shown by Damanik and Lenz |DL99a| that, given any w £ ^a,b and fc > 1, 
the restriction of to {rt > 1} begins with a square 

Wl . . . W2Ffe . . . = CJl . . . LOPk^l ■ ■ ■ W_Ffc ■ • ■ 

such that uJi . . . uip^ is a cyclic permutation of S'' (a) . By cyclic invariance of the 
trace, it follows that TTMi^{Fk,E) — 2xk{E) for every E. 

The Cayley-Hamilton Theorem, applied to M^{Fk,E)^ says that 

M^{Fk,Ef - (TrAf„(Ffe, S)) Af^(Ffc, E) + I = Q, 

which, by the observations above, translates to 

M^(2Ffc, E) - 2xt,MjFk, E) + I = 0. 



If G Sa,6 and M is a solution of the difference equation (47), it therefore follows 
that 

U{2Fk + 1) - 2xkU{Fk + 1) + U{1) - 0. 

If u does not vanish identically, this shows that u„ 7^ as n — > oo since the Xk 's are 
bounded above and the w„'s are bounded below away from zero. In particular, if 



E G then no non-trivial solution of (47 1 is square-summable at +00 and hence 



E is not an eigenvalue. It follows that the point spectrum of is empty. □ 
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